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and  compared  with  the  thermometer  record  for  those  years  to  obtain 
phenomenological  (numerical)  temperature  coefficients  for  the  ratios 
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SECTION  1.  RESEARCH  PROGRAM 


Mechanical  thermometers  were  invented  about  1680  A.D.  The  records  of 
temperature  for  the  next  50  years  or  so  are  difficult  to  correlate 
with  each  other,  however,  because  the  calibrations  differed  from  ther¬ 
mometer  to  thermometer.  Nevertheless,  heroic  efforts  have  been  made 
to  do  so,  and  thus  the  record  of  the  long,  slow  increase  in  air  tempera¬ 
ture  since  the  Little  Ice  Age  has  been  unraveled.  If  one  wishes  to 
study  climate  changes  before  this  time  it  is  necessary  to  invent  and 
calibrate  other  kinds  of  thermometers. 

Some  25  years  ago  Harold  Urey  made  the  first  attempt  by  measuring  the 

isotope  effect  in  the  carbonate  of  shells,  that  is,  the  variation  of  the 
13  12  18  16 

ratio  of  C  /C  or  0  /0  in  shells  grown  in  water  of  different  tempera¬ 

tures.  (The  values  of  the  ratio  are  larger  at  higher  temperature.)  However, 
in  order  to  evaluate  the  water  temperature  in  which  the  shell  grew,  it  is 
necessary  ic  know  the  isotope  ratios  in  the  water  and  in  the  carbonate- 
bicarbonate  ion  where  the  shell  grew.  Since  in  general  these  are  not  known 
for  old  shells,  the  isotope  thermometer  in  shells  has  not  been  very  useful. 
Furthermore,  this  method  provides  only  one  thermometer,  because  carbon  and 
oxygen  are  bound  through  the  same  force  constant  k,  and  therefore  vibrate 
against  each  other  with  the  same  frequency,  u.  Since  the  isotope  ratio 
depends  only  on  hw/kT,  the  temperature  effects  on  carbon  and  on  oxygen  are 
proportional  to  each  other  and  hence  give  redundant  information. 

In  the  present  research,  the  isotope  thermometers  in  bio-organic 
material  are  being  calibrated.  Here  there  are  five  independent 
thermometers,  namely  hydrogen,  carbon,  nitrogen,  oxygen  and  sulfur. 


Any  given  one  of  these  temperature  coefficients  should  have  the  same 
numerical  value  in  all  bio-organic  materials  because  its  bonds  have 
about  the  same  strengths  in  all  bio-organic  molecules.  Ir  parallel 
changes  are  observed  in  two  or  more  of  these  thermometers  in  stored 
bio-organic  material  from  past  times,  it  becomes  increasingly  likely 
that  a  climate  change  i.3  the  cause.  Accordingly,  this  research  program 
became  oriented  toward  measuring  temperature  coefficients  of  bio-organic 
material  both  in  general  -.ituations  where  atmospheric  C02  (which  has 
the  same  isotopic  composition  world  wide)  and  rain  water  are  involved, 
and  also  in  special  fresh-water  environments.  The  ramifications  and 
applications  appear  to  be  large. 

At  present,  the  climate  of  the  northern  hemisphere  appears  to  be  in  a 
long-term  decline  that  began  in  1935.  Accompanying  the  decline  in 
temperature  have  been  50  year  droughts  in  Japan,  North  Africa,  India, 

Scuth  America  and  the  Mediterranean.  Indicative  of  this  trend  is  a 
statement,  by  the  head  of  the  Soviet  weather  service.  Dr.  Yevgeny  K. 
Federov  quoted  in  the  New  York  Times  that  temperatures  in  1972  in 
central  Russia  w’.re  the  lowest  in  several  hundred  years.  If  the 
temperature  decline  and  the  long-term  droughts  persist  the  effect  on 
world  food  supply  could  become  critical.  In  illustration  of  the  recent 
rapid  decrease  in  world  food  reserves,  is  an  estimate  of  Lester  R.  Brown, 
Overseas  Development  Council,  New  York  Times,  January  26,  1974,  of 
reserves  for  95  days  in  1961,  foj  51  days  in  1971,  and  for  29  days  in 
January  1974.  This  estimate  presumably  takes  into  account  the  increase 
in  world  population  of  about  30%  during  the  interval. 

The  food  shortage  is  likely  to  be  enhanced  by  the  fuel  shortage  because 
this  has  produced  a  shortage  of  fertilizer  on  which  the  new  "Green  Revolu¬ 
tion"  high  yield  grains  are  dependent,  being  unable  to  thrive  on  old- 
fashioned  animal  manures. 
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The  present  research  program  may  succeed  in  developing  a  method  to 
evaluate  temperature  trends  in  the  historic  past,  say  for  the  last 
10,000  years,  and  from  the  record  to  determine  the  pattern  of  climate 
changes,  from  which  we  may  be  able  to  better  predict  the  climate  changes 
to  be  expected  in  the  near  term  and  their  potential  impact  on  food 
supply. 


SECTION  2.  MAJOR  ACCOMPLISHMENTS 


he  following  references  and  abstracts  describe  the  studies  performed 
under  this  contract. 

"Measurement  of  0^/0^  by  Neutron  Activation,"  by  R.  G.  Jackson, 

L.  M.  Libby,  and  H.  R.  Lukens,  J.  Geophys.  Res.,  78,  p.  7145  (1973). 

(See  Appendix  A. ) 

The  commercially  available  services  for  isotope  measurement  of  an 

accuracy  useful  in  evaluation  of  past  climates  from  bio-organic  material 

are  exceedingly  limited.  In  seeking  to  alleviate  this  shortcoming,  and 

in  accordance  with  one  of  the  tasks  of  this  contract,  we  have  proven 
18  16 

that  the  ratio  0  /0  can  be  measured  in  water  using  neutron  activation 

4 

to  an  init:'al  accuracy  of  a  few  parts  in  10  ,  with  good  prospects  of 
improved  accuracy.  By  a  small  increase  in  detector  sensitivity  the 
method  promises  to  become  of  commercial  importance. 

"Isotope  Ratios  in  a  Bristlecone  Pine,"  L.  M.  Libby  and  L.  J.  Pandolfi. 

18  16  13  17 

We  measured  the  ratios  0  /0  and  C  /C  '  in  a  piece  of  bristlecone 
pine  for  the  years  500  B.C.  to  200  A.D.,  with  the  intention  of  comparing 
them  with  those  of  a  recent  pine  grown  at  the  same  site  (Almagre 
Mountain,  Colorado  Springs,  Colorado,  11,000  ft  altitude).  Then  we 
learned  by  study  of  the  recent  pine  that  each  tree  ring  remained  in  the 
sapwood  for  as  long  as  90  years,  permitting  isotope  exchange  between 
cellulose  and  sap  over  this  period  and  possibly  smearing  any  tempera¬ 
ture  record  over  this  period.  The  temperature  record  made  at  the  site 
from  mercury  thermometer  readings  extends  over  only  the  last  90  years. 
Thus  we  could  not  calibrate  the  observed  isotope  changes  using  this 
record.  However,  we  can  say  that  for  both  018/016  and  C13/C12  ratios 
the  difference  between  the  old  pine  and  the  modern  pine  indicated 
that  the  climate  was  warmer  2000  years  ago  than  now. 


The  measurements  are  shown  in  Table  1.  The  age  of  the  old  tree  was 
determined  by  radio  carbon  dating,  by  Dr.  Rainer  Berger,  UC1A,  to  a 
precision  of  +60  years.  The  modern  tree  rings  were  dated  by  Dr.  Valmont 
La  Marche,  Tree  Ring  Laboratory,  University  of  Arizona,  who  cored  the 
tree  and  courted  the  rings. 


Table  1.  Average  Ratios  in  Bristlecone  Tree  Rings, 
1896-19)51  A.D.  and  550  B.C.-75  A.D. 


(C13/C12)* 

(018/016)* 

Time  Span 

Modern 

Old 

.013062 

.013240 

.0024924 

.0025400 

1896-1951  A.D. 

550  B.C.-75  A.D. 

*With  errors  of  about  0.5  parts  per  thousand 


As  shown  in  Table  1,  both  isotope  ratios  are  larger  in  the  old  tree  than 
in  the  modern  tree.  Based  on  these  data,  thermodynamic  considerations 
indicate  that  the  old  tree  grew  in  a  wanner  climate.  But  since  there 
is  available  no  experimental  determination  of  the  temperature  coefficient 
for  these  ratios  in  Bristlecone  pines,  we  can't  say  how  much  warmer. 
However,  if  one  assumes  that  the  Bristlecone  temperature  coefficients 
may  be  equal  to  these  we  have  evaluated  for  the  E  ropean  Spessart  Oak, 
of  4.4  +0.3  ppt/°0  for  oxygen  and  1.8  +0.3  ppt/°C  for  carbon,  the 
dffference  in  the  oxygen  ratio  would  indicate  a  climate  colder  today 
by  4.3  +0. 3°C  and  the  carbon  difference  would  indicate  7.6  +1.4°C  colder, 
for  an  average  climate  change  of  6.0  +1.4°C.  Nevertheless,  there  is  no 
reason  to  believe  in  the  validity  of  this  number,  both  because  tempera¬ 
ture  coefficients  for  pine?  may  differ  from  those  for  oaks  and  also 
because  these  numerical  temperatures  coefficients  (4.4  and  1.9  ppt  per 
°C)  were  evaluated  from  a  smaller  temperature  change,  about  1.8  degrees. 
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If  the  temperature  coefficients  are  not  independent  o'  temperature, 

there  may  be  a  substantial  error  in  using  them  to  evaluate  a  temperature 
change  of  several  degrees. 

It  is  clear  that  many  new  questions  have  arisen  in  the  course  of  this 

research  project;  it  is  important  to  evaluate  the  answers  in  order  that 

this  potential  new  field  of  historic  climate  determination  may  become 
viable. 


"Temperature  Coefficients  of  Stable  Isotopes  Measured  in  Hot  Spring 
Algae/'  by  L.  M.  Libby  and  L.  J.  Pandolfi,  is  appended  herewith, 
together  with  the  first  page  of  proof  for  its  publication,  (see 
Appendix  B.) 


The  algae  in  which  isotope  ratios  were  measured  grew  at  various  places 
of  successive  lower  temperatures  in  the  downflow  of  the  famous  hot 
springs  of  Rotorua,  New  Zealand.  In  evaluating  the  temperature  coef¬ 
ficients  of  the  isotope  ratios  of  these  algae,  we  have  assumed  that 
the  isotopic  ratio  of  the  spring  water  remained  unchanged  during  the 
flow  down  hill.  To  check  this  point  I  have  written  to  Dr.  Athol  Rafter, 
Director  of  the  Institute  of  Nuclear  Sciences,  and  to  Mr.  Thomas  Grant-’ 
Taylor,  Geological  Survey,  both  in  the  Department  of  Scientific  and 
Industrial  Research,  Lower  Hutt,  New  Zealand,  asking  if  they  will  have 
the  isotope  ratios  measured  in  water  samples  obtained  at  60°  and  at 
30°  centi8^e  on  the  downhill  slope  or  alternatively  send  me  the  samples 

"Experimental  Determination  of  Temperature  Coefficients  of  Isotope 
Ratio-  in  Elodea,"  by  L.  M.  Libby  and  L.  J.  Pandolfi. 


The  angiosperm  Elodea  was  grown  in  an  outdoor  pond,  temperature  controll 
at  10°C  and  20°C.  The  crop  of  new  tips  was  harvested  after  about  2  mont 
of  growth  at  each  temperature  and  the  C13/C12  ratios  measured  in  the 


6 


material.  To  effect  these  measurements,  the  entire  tips  were  dried 
in  vacuo,  and  then  subjected  to  chemical  treatments  for  evolution  of 
gases  measureable  in  a  mass  spectrometer.  (See  Appendix  C.)  The  data 
are  shown  in  Table  II  as  well  as  the  experimentally  evaluated  tempera¬ 
ture  coefficient. 


Table  2.  Isotope  Ratio  Measured  in  Pond  Grown  Elodea 


Water 

Temperature 

c13/c12 

Experimental 

Temperature 

Coefficient 

4°C 

12992  +43 

10°C 

13009  +38 

(0.26  +  0.13)  ppt/°C 

20°C 

_ 

13053  +55 

During  the  course  of  this  experiment,  because  the  pond  was  outdoors 
and  quite  small,  9  cu  ft,  it  was  subject  to  rapid  evaporation,  so  that 
make-up  water  had  to  be  added  almost  every  day,  from  the  Los  Angeles 
water  supply.  The  source  of  the  Los  Angeles  water  supply  is  in  part 
the  melting  snows  of  Mammoth,  whose  isotope  ratios  vary  with  the  time 
of  year  at  which  the  snow  is  deposited  because  the  air  temperature 
of  the  snow-clouds  is  variable.  Thus  the  isotope  ratio  of  the  make-up 
water  changed  during  the  experiment  as  it  was  produced  from  the  melt¬ 
ing  snows.  On  account  of  this  varying  and  uncontrollable  effect,  we 
could  not  correct  the  measured  oxygen  and  hydrogen  isotope  ratios  in  a 
reliable  way  and  so  can  not  report  our  measurements.  To  evaluate  their 
coefficients  will  require  growing  plants  under  conditions  where  the 
make-up  water  comes  from  a  water  tank  so  enclosed  that  evaporation  can¬ 
not  occur,  and  of  sufficiently  large  volume  as  to  last  for  the  duration 
of  the  experiment. 


"Calibration  of  Isotope  Thermometers  in  a  European  Oak,"  L.  M.  Libby 
and  L.  J.  Pandolfi,  presented  at  the  NSF  International  CLIMAP  Meeting 
on  Climatic  Parameters  Since  17,000  B.C.,  May  17-23,  1973,  University 
of  East  Anglia,  Norwich,  England,  in  publication  in  the  Proceedings; 
and  presented  at  the  International  Colloquim,  C.N.R.S.,  Centre  des 
Faibles  Radioactivite,  Oef-sur-Yvette,  France,  "Quantitative  Methods 
of  Study  of  Climate  in  the  Pleistocene,"  June  5-8,  1973,  and  published 
in  the  Proceedings  (first  page  of  proof  appended  herewith),  further 
reviewed  by  Prof.  Hans  Suess,  UCSD  (See  Appendix  C). 

We  measured  the  isotope  ratios  of  carbon,  hydrogen,  and  oxygen  in 
rings  of  a  German  oak  for  the  years  1710  to  1950  and  compared  them 
with  measured  air  temperatures  of  England,  Basil  and  Geneva  for  the 
same  time  span.  In  this  way,  we  obtained  numerical  evaluations  of  the 
temperature  coefficients  of  all  three  isotope  ratios. 

The  isotope  ratios  have  been  measured  in  an  oak  obtained  from  Dr.  B. 

Becker,  Landwirtschaf tliche  Hochschule,  Universitat  Hohenheim,  Stuttgart. 

It  was  grown  in  the  Spessart  mountains  in  1514-1850  A.D. :  its  rings 

have  been  cross-dated  with  the  Huber  master  tree  ring  chronology  by 

Becker.  The  measured  ratios  are  shown  versus  years  in  Figure  1.  The 

13  12 

temperature  coefficient  for  C  /C  ,  as  determined  by  the  warm  period 

around  1730,  is  comparable  with  that  measured  for  the  oak  from  Munich. 

18  16 

However,  the  coefficients  for  D/H  and  0/0  seem  to  be  only  half  as 
large  as  for  the  Munich  oak,  suggesting  that  this  second  oak  was 
partly  nourished  by  old  ground  water  which  was  insensitive  to  isotope 
changes  in  yearly  precipitation.  However,  Becker  has  no  information 
on  the  source  of  water  for  the  second  oak.  The  records  for  all  three 
indicate  that  the  years  1540-1700  were  cold,  in  agreement  with  what 
is  known  about  the  Little  Ice  Age. 
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I  sent  the  data  for  this  second  oak  (Figure  1)  to  Professor  H.  H.  Lamb, 
Director  of  the  School  of  Environmental  Sciences,  University  of  East 
Anglia,  Norwich,  Norfolk,  England,  for  his  comments.  He  replied  as 
follows : 


13 

"The  C  temperature  curve  on  the  Spessart  oak  looks 
to  me  as  if  it  correlates  well  with  the  variations 
of  Baltic  ice  (dates  of  opening  of  the  port  of  Riga). 
All  your  curves  show  the  interesting  warm  intervals 
during  the  cold  ("Little  Ice  Age")  time,  viz.  mid 
1600s,  also  1570s,  1730s,  1790s  with  which  I  (i.e.  my 
various  data)  agree." 
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Figure  1.  Isotope  Ratios  Measured  in  the  Second  Spessart  Oak 
Plotted  Versus  Dates  Evaluated  From  Official  Tree 
Ring  Analysis. 
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2,1  CLIMATIC  INFLUENCES  ON  RADIOCARBON  DATES 

This  portion  of  the  research  program  has  the  objective  to  account  for 
the  complete  carbon  inventory  and  its  perturbations.  The  following 
reports  have  been  prepared  to  date: 

a)  "Vulcanism  and  Radiocarbon  Dates,"  by  L.  M.  Libby  and  W.  F. 
Libby ,  Proceedings  of  the  9th  International  Radiocarbon  Dating 
Conference.  October  1972,  Lower  Hutt,  New  Zealand,  ed.  T.  A. 
Rafter  and  T.  L.  Grant-Taylor ,  page  A72  (See  Appendix  D). 

We  compared  the  amount  of  CO2  emitted  from  volcanoes  with  the 
amount  of  CO  2  being  injected  into  the  atmosphere  by  burning 
of  fossil  fuels  we  are  able  to  conclude  that  volcanoes  cannot 
cause  important  perturbations  unless  they  Increase  in  frequency 
by  a  factor  100  or  more. 

b)  "Production  of  Radiocarbon  in  Tree  Rings  by  Lightningbolts," 
by  L.  M.  Libby  and  H.  R.  Lukens,  J.  Geophvs.  Res..  78,  p.  5902 
(1973)  (See  Appendix  E). 

Consideration  of  neutro;.  production  in  electric  discharge 
through  threads  in  the  laboratory  gives  an  experimental  basis 
for  estimating  neutron  production  in  lighting  discharges. 
Radio-carton  production  in  a  tree  struck  by  lightning  may  be 
measureatle,  of  order  1%,  of  the  radio  carbon  produced  by 
cosmic  ray  neutrons. 

c)  "Globally  Stored  Organic  Carbon  and  Radiocarbon  Dates,"  by 

L.  M.  Libby,  J.  Geophys.  Res.,  78,  p.  7667  (1973)  (See  Appendix 
F). 
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Using  recent  data  on  the  rate  of  addition  of  C02  to  the  atmos¬ 
phere  by  bacteria  it  is  found  that  variations  of  climate  acting 
on  the  biosphere  can  cause  fluctuations  of  up  to  +100  years 
correction  to  radio  carbon  dates. 
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SECTION  3.  CONCLUSIONS 


It  is  perhaps  the  basic  nature  of  fundamental  research  that  new  problems 
arise  as  a  consequence  of  trying  to  answer  an  initial  problem.  In  the 
present  research  a  large  number  of  new  questions  have  appeared,  such  as 

1)  are  the  isotope  temperature  coefficients  alike  in  all  European  oaks? 

2)  are  other  kinds  of  trees  similar  to  oaks  in  that  respect,  in  parti¬ 
cular  the  Bristlecone  pines  for  which  there  exist  chronological  sequences 

tree  rings  for  the  last  8000  years?  and  3)  are  the  temperature 
coefficients  for  isotope  fractionation  in  manufacture  of  cellulose 
and  for  formation  of  rainfall  separately  evaluable?  The  answer  to  this 
third  question  is  basic  to  understanding  Isotope  thermometry  in  trees. 

To  evaluate  these  two  effects  separately  we  must  grow  plants  under 
conditions  where  both  temperature  and  the  isotope  ratios  of  the  water 
supply  are  controlled.  Under  these  conditions,  the  measured  temperature 
coefficients  of  the  Isotope  ratios  will  be  those  characterizing  manu¬ 
facture.  The  temperature  coefficients  measured  for  tree  rings  in 
natural  environments  are  equal  to  the  sum  of  the  coefficient  for  chemi¬ 
cal  manufacture  of  cellulose  plus  the  coefficient  for  formation  of 
precipitation  in  the  troposphere  above  the  tree. 

Numerical  results  to  date  are  summarized  in  Table  III. 


Consultants  to  this  project  are  Prof.  Hans  Suess  and  Harmon  Craig 
(UCSD) ,  V.  C.  LaMarche  and  B.  Bannister  (University  of  Arizona),  R. 
Berger,  L.  Pandolfi  and  W.  F.  Libby  (UCLA),  H.  G.  Jackson  and  J.  Enstrom 
(UCB) ,  and  S.  Korff  (NYU). 
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Table  III 


Summary  of  Temperature  Coefficients  of  Isotope 
Fractionation  in  Organic  Materials 


A. 


43 


Temperature  Coefficients  for  Biological  Materials  (ppt/°C) 
Measured  for  Temperature  Ranges  of  a  Few  Degrees  C. 


Wood 


Atmospheric  ^ 

1.8  +0.3 


Marine 

Phytoplankton  0.35 

Plankton  0.5 


Hot  Springs 

Bluegreen  Algae  -0.62  +0.12 

Los  Angeles  Water  Supply 
Eu«lena  3.9  +0.6 


Atmospheric  Moisture 
Wood  4.4  +0.3 

Hot  Springs 

Bluegreen  Algae  2.2  +1.0 


Libby  and  Pandolfi 

Degens  et  al. 

Eadie 

Libby  and  Pandolfi 
Libby  and  Pandolfi 

Libby  and  Pandolfi 
Libby  and  Pandolfi 


C.  D 


Wood 

Peat 


Atmospheric  Moisture 
95  +19 

Ground  Water 

3 


Libby  and  Pandolfi 
Schiegl 
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Table  III  (Continued) 


Measured  Isotope  Ratios  in  Algae,  Expressed 
in  Parts  per  Thousand,  Vs.  Water  Temperature  in 
Degrees  Centigrade 


6  * 

T  °13 

6  * 

°18 

degrees 

parts  per 

parts  per 

Centigrade 

thousand 

thousand 

u> 

o 

o 

+13.7 

— 

46 

-  2.6 

-23.3 

51 

+  1.4 

+17.0 

55 

-  1.2 

+15.1 

56 

-  3.3 

+17.7 

64 

-  9.3 

+22.4 

*6^  and  5^g  are  the  differences  from 
50°C  by  least  squares  analysis^. 

613  =  (31.1  -  [0.62  40. 12]T)  ppt 

6lg  *  {“HO  +  [2.2  +1.0]t}  ppt 

the  values  pre 

Temperature  Coefficients  for  Biological  Materials  Computed  for 

Cellulose  ppt/°C 


a18(02) 

0.96 

a18(C02) 

1.14 

a18(H20) 

0.92 

13 

0.36 

a 

a^HCOD) 

0.4 

a(DCOH) 

2.0 

**See  Appendix  G. 
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APPENDIX  A 

MEASUREMENT  OF  018/016  RATIO  USING  A  FAST  NEUTRON  REACTOR 
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VOL.  78,  NO.  30 


JOURNAL  OF  GEOPHYSICAL  RESEARCH 


OCTOBER  20.  1972 


Measurement  of  tH0/'*0  Ratio  Using  a  Fast  Neutron  Reactor 

H.  G.  Jackson 

Lawrence  Livermore  Laboratory,  Livermore,  California  91,560 


L.  M.  Libby 

School  of  Engineering  and  Applied  Sciences,  University  of  California  at  Los  Angeles 
Los  Angeles,  California  90021, 

H.  R.  Lukens 

Gulf  Radiation  Technology,  La  Jolla,  California  92037 

A  method  has  been  devised  and  proven  to  measure  the  stable  isotope  ratio  ,'0/'*0  for 
oxygen  by  neutron  activation  in  a  fast  neutron  reactor.  The  measurement  is  made  non- 
destrur lively  on  a  few  grams  of  water  with  an  neruraey  of  a  few  parts  in  10,000.  A'  present, 
measurement  to  this  accuracy  requires  about  40  min  with  the  reactor  running  at  a  power  of 
250  kw,  but  by  increasing  detector  sensitivity  the  method  promises  to  become  a  routine  matter 
accomplished  in  less  time. 


If  rapid  nondestructive,  ways  of  measuring 
ratios  of  stable  isoto|M'.s  were  available,  the  use 
of  stable  isotopes  in  research  would  be  greatly 
aided.  Wc  have  shown  that  for  the  stable 
isotopes  of  oxygen  the  ratio  " 0/"0  ran  be 
measured  nondot ructivelv  in  small  amounts  of 
water  by  neutron  activation,  he  activities  inea- 
sured  are  {a)  "0(«.  ;»)K'X;  threshold  neutron 
energy  is  10.25  Mcv.  (<7„>(«.  ;>)  =  40  mb.  tH 
(mean  life)  =  10, '1  see,  and  /?,  =  0,14  Mcv 
(OSff )  and  ( b )  “0(/i,  y ) '"O ;  formed  by  thermal 
neutron  capture,  <r,dy,  *t)  =  0.22  mb,  t„  (mean 
life)  =  41.8  sec,  A'„  =  1.37  Mcv  (60%),  and 
=  200  kev  (100%), 

It  is  necessary  to  irradiate  the  water  in  a 
neutron  reactor  that  has  a  high  fractional  flux 
of  fast  neutrons  (relative  to  that  available  in 
the  graphite-moderated  reactors)  in  ortler  for 
the  10-sec  activity  to  have  .sufficient  intensity. 
Its  gamma  ray  has  a  high  energy.  0.14  Mcv, 
and  so  is  relatively  easy  to  measure,  because 
the  background  activities  from  impurities  in 
water  and  from  the  Compton  radiation  are  small 
at  this  energy.  The  42-see  activity  has  two 
gamma  rays,  namely,  1.37  and  0.2  Mcv.  That 
at  1.37  Mcv  would  be  the  more  desirable  to 
measure  except  that  it  coin  rides  with  a  gamma 

Copyright  ©  1973  by  ihe  American  Gi  m  ay  si  cat  Union. 


ray  from  neutron -activated  sodium  (which  is 
an  abundant  impurity)  so  closely  that  it  cannot 
be  identified  separately.  Therefore  this  activity 
was  measured  by  counting  the  lower-energy 
gamma  radiation  at  200  kev.  Here  die  Compton 
background  i-  large  and  is  die  factor  limiting 
the  accuracy  of  each  determination. 

Thr  10-  and  42-sec  activities  are  measured 
at  the  same  time  in  the  same  detector,  being 
separated  by  pulse  height  analysis,  and  recorded 
separately.  Water  samples  in  which  it  is  desired 
to  measure  the  ratio  "0/"0  are  compared  with 
a  standard  water  sample,  and  the  ratio  relative 
to  the  standard  is  expressed  as  a  deviation  8” 
in  parts  per  thousand  (%„)  according  to  the 
definition 

<n?<)  =  l(rt. . /«., -  lj  X  1000  (1) 

The  desired  accuracy  is  obtained  by  irradiating 
and  counting  the  induced  activities  of  the  water 
sample  several  times  until  some  hundred  million 
gannna  rays  have  been  counted. 

Kxr Kill  MENTAL 

The  Triga  Mark  I  reactor  in  steady  state 
njieration  at  250  kw  has  a  thermal  neutron  flux 
of  2.8  •  lO'Ycm'  .-ee  and  a  fast  (fission  spec-, 
trnm)  flux  of  3.5  •  I0ls/ctn3  sec  at  the  irradia¬ 
tion  position.  Water  samples  of  approximately 
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5  grams  each  in  plastic  vials  (purchased  from 
Olympic  Plastics  Company,  Los  Angeles,  Cali¬ 
fornia)  containing  no  intrinsic  oxygen  were 
shot  into  the  reactor  in  the  fast  pneumatic 
transfer  system  and  there  irradiated  for  30  see, 
after  which  they  were  shot  out  into  a  repro¬ 
ducible  site  in  the  counting  facility.  After  an 
interval  of  10  sec  the  induced  10-  and  42-sec 
activities  were  counted  during  the  next  subse¬ 
quent  interval  of  30  see. 

Reactor  power  is  maintained  at  a  present 
level  hv  servo-regulated  control  rod  drives.  The 
insertion  of  5  grams  of  water  in  plastic  vials 
into  the  outermost  fuel  element  ring  (sample 
irradiation  jiosition)  is  equivalent  to  adding  a 
slight  amount  of  reactivity  to  the  core,  which  is 
immediately  compensated  for  by  automatic  con¬ 
trol  rod  adjustment  without  overshoot.  The 
reactor  |iower  is  maintained  to  ±0.13%  of  the 
preset  level  at  a  confidence  level  of  95%.  The 
tem|>erature  at  the  sample  position  ranges 
from  20°  to  25#C  during  the  year  but  is  held 
ron>f.mt  during  any  given  hour  of  reactor 
operation  bv  forced  circulation  of  pool  water. 
The  technology  and  use  of  the  reactor  in  simi¬ 
lar  measurements  were  previously  described 
[Bramblett  et  al.,  1973], 

The  count  in  the  200-fcev  channel  was  at  a 
rate  of  0.2  megacycles,  and  this  is  too  fast  for 
standard  pulse  height  selection  equipment. 
Therefore  a  simple  pulse  height  analyzer  was 
huilt  (Figure  1)  with  integrated  circuitry.  The 
system  is  limited  to  a  maximum  count  rate  of 


300  megacycles  by  the  sealers.  The  sensor  was 
a  Nal  crystal  %  inch  high  and  2  inches  in  di¬ 
ameter  glued  with  a  high-viscosity  fluorocarbon 
(Dow  Corning  20-057)  to  a  photomultiplier, 
RCA-6810  A. 

Tho  counts  observed  in  5  grams  of  normal 
wafer  (Colorado  River  water  from  the.  public 
water  system  of  La  Jolla,  California,  for  which 
“0/‘0  was  measured  in  a  mass  spectrograph) 
in  a  counting  interval  of  30  sec  for  the  200-kev 
radiation  (channel  2)  and  the  6.1-Mev  radia¬ 
tion  (channel  4)  are  listed  in  Table  1.  Counting 
ratios  measured  at  two  different  power  levels  of 
the  reactor,  namely,  at  190  and  250  kev,  are 
strictly  comparable,  showing  that  counting  losses 
were  negligible.  Varying  amounts  of  water  en¬ 
riched  in  "O  were  added  to  normal  water  to 
prepare  a  series  of  samples  in  which  "O  was 
progressively  concentrated.  The  measured  ratios 
of  counts  in  the  two  channels  recording  the  200 
kev  and  the  6.0  Mev  radiation,  respectively,  are 
listed  in  Table  2  for  increasing  concentrations 
and  are  shown  to  be  linearly  dependent  on  the 
"O  concentration.  As  is  shown  in  Table  1,  after 
24  runs,  each  of  70  see,  the  ,,0/uO  ratio  in  a 
given  water  sample  was  measured  to  an  overall 
error  of  ±0.06%,  after  24  determinations,  each 
of  standard  deviation  0.29%.  Xo  external  per¬ 
turbations  were  allowed  to  influence  the  reactor 
during  the  measurements.  The  overall  standard 
deviation  (0.29%/ (24)-  =  0.06%,)  testifies  to 
the  constancy  of  the  ratio  of  fast  neutron  flux 
[Dixon,  1971]. 


Fig.  1.  Block  diagram  of  pulse  height  selector. 
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TA8LE  1.  Observed  Decays  Proportional  to  Concentrations  of  ‘80  and  of  ,60  and  Their  Ratios 
for  24  Experiments  on  Faucet  Water  from  the  Public  Watev  Supply  of  La  Jolla,  California 


Sample 

Number 

Net  Weight, 
grams 

Net  Counts 

x,  Count  Ratio 
Channel  2/Channel  4 

Channel  2 

Channel  4 

Sept. 

i  7,  190-h)  Rtictor  Pcrjer 

1 

4.9640 

3,793,063 

531,946 

7.130S 

2 

4. 9643 

3,805,560 

536,306 

7.0960 

3 

4.9639 

3,872,147 

546,752 

7.0821 

4 

4.9681 

3,756,711 

.'.31,097 

7.073S 

S 

4.9648 

3,600,182 

507,927 

7.0880 

6 

4.9694 

3,734,811 

528,492 

7.0671 

7 

4.9661 

3,675,843 

515,083 

7.1364 

8 

4.9672 

3,771,143 

531,287 

7.0981 

9 

4.9718 

3,708,926 

523,682 

7.0824 

10 

4.9666 

3,788,234 

SV  952 

7.0814 

a 

4.9695 

3,814,125 

538,040 

7.0890 

12 

4.9673 

3,758,823 

530,097 

7.0908 

Sept. 

12,  250~kw  Reactor  Power 

1 

4.9640 

4,123,664 

580,767 

7,1004 

2 

4.9643 

4.0S7.737 

571,263 

7,1031 

3 

4.9639 

4,309,732 

608,763 

7.0795 

4 

4.9681 

3,975,696 

560,343 

7.0951 

5 

4.9648 

4 ,295,982 

605,418 

7.0959 

6 

4.9694 

4,210,863 

597,017 

7.0531 

7 

4.9661 

4,230,811 

598,442 

7.0697 

8 

4.9672 

3,884,063 

550,786 

7.0518 

9 

4.9718 

3,884,013 

549,599 

7.0670 

10 

4.9666 

3.874,433 

548,264 

7.0667 

11 

4.9695 

4,331,802 

612,979 

7 . 0668 

12 

4.9673 

4,074,085 

574,198 

7.0952 

JO-sec  irradiation,  10-sec  delay,  JO-sec  count. 
<*>  •  7. 0858;  o(4)  -  0.29;  U<1>(\)  •  0.0S9. 


ob,*n'e<1  0'C*V5  Proportion*!  to  Concentre! ionj 
of  0  and  of  “0  in  Veter  Sample]  Se.uentiellp  Enriched 
In  ‘*0  vertui  >»0/“o  Coapoted  fro.  Mess 
Spectroscopic  Keesurements 


Computed  X 

Jl0,  wtt  of 

Tot  el  Oxygen* 

Meesured  1 

Channel  Retio, 
Channel  2/Channel  4 

t  Predicted  from 
Regression  line* 

0.22044 

0.25175 

0.24594 

0.24754 

0,27510 

.27719 

0.50902 

0.507SS 

0.54506 

0.5S92S 

7.0929 

7  1687 
.  S79 

7.2592 

7.4089 

7  4129 

7.6206 

7.S950 

7.7980 

7.8860 

7.1050 

7.1692 

7.2497 

7.2587 

7.4190 

7.4270 

7.6974 

7.5991 

7.8005 

7.8924 

Pearson  product .Moment  correlation  cooffieiont  rr»  • 
0.9-J9JS;  "prediction  *  tO-l.’M  roletlv.; 

"■ten  prediction  *  lO.Oait  relative. 

All  but  first  sample  prcpercd  by  addition  of  enriched 
water.  All  but  first  semple  meesured  once.  First  semplt 
Le  Jo‘le  tep  water,  wes  measured  24  times  (Teble  If. 

♦Regression  line  Is  obtained  from  observables:  /  • 

S. 6721426(f)  *  S.IS4646.  The  lest  number  is  the  ratio  of 
backgrounds  in  channel  2  to  channel  4  in  tho  ebsence  of 


Tho  enriched  water  was  measured  by  mass 
sjicctrograph  [Ritterberg  and  Pontecorvo, 
1955], 

Tho  pneumatic  transfer  apparatus,  Triga 
Mark  I  reactor,  and  timing  system  have  been 
used  as  unrl.ingcd  components  in  a  “U  mea¬ 
surement  procedure  that  delivered  0.37%  pre¬ 
cision  in  individual  measurements  [Bramblctt 
rt  a!.,  1973].  Hence  the  variation  in  individual 
channels,  a  ~  i‘2%  on  September  7  and  a  S 
■1%  on  September  12,  is  believed  to  be  controlled 
by  the  "'"inting  equipment.  Fortunately,  the 
couni  ratio  proved  to  be  largely  independent 
of  the  variance  of  the  individual  channel.  How¬ 
ever,  it  is  likely  that  improvements  in  counting 
equipment-  stability  anil  counting  ellirieneies 
will  improve  both  the  speed  and  the  precision 
of  the  isotojK*  ratio  measurement. 

Conclusion 

A  nondestructive  method  of  measuring  the 
stable  isotope  ratio  “0/wO  in  water  has  been 


A-4 


Jackson  et  al, 

sliown  to  be  possible,  rapid,  and  feasible  as  a 
commercial  service,  to  be  presently  accurate  to 
about  4  parts  in  10,000  at  a  reasonable  ex¬ 
penditure  of  measuring  time  (about  40  min  at 
present),  and  to  promise  fiitur'  measurements 
at  least  to  this  accuracy  in  s  ior'er  times,  whtn 
the  detector  sensitivity  has  been  increased. 
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MEASURED  TEMPERATURE  COEFFICIENTS 
OF  STABLE  ISOTOPES  IN  HOT  SPRING  ALGAE 
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Abstract 

The  temperature  coefficient  of  three  stable  isotopes  in  a  blue-green  algae,  growing  in  alkaline  hoi  springs  at  temperatures 
of  30  C-68  C.  have  been  measured.  That  for  oxygen  is  found  to  be  +(2.2*  1.0)  ppt/  'C  which  is  larger  than  the  value 
computed  for  equilibrium  as  is  to  be  expected.  That  for  carbon  is  found  to  be  (—0.62*0.12)  ppt/"C;  the  negative 
value  is  attributed  to  methane  formation  by  copious  bacteria  so  that  the  system  CH,  —  HCO.,  which  is  known  to  have 
a  negative  temperature  coefficient,  dominates.  That  for  hydrogen  is  found  to  be  (3.3  *  0.7)  ppt/  'C,  about  twice  that 
computcJ  lor  equilibrium  lormation  of  cellulose.  (The  quoted  errors  represent  one  standard  deviation). 


Rlsu.mC 

Lcs 'Coefficients  do  temperature  do  trois  isotopes  stables  pour  dcs  algucs  blcu-vcrt  sc  developpant  dans  dcs  sources 
chaudes  a  des  temperatures  dc  3tl  "C  a  68  "C  ont  tie  mesurcs.  On  a  trouve  que  Ic  rapport  pour  l’oxygenc  est  2,2  *  1.0  ..  "C 

ce  qtu  cm  plus  eleve  que  la  v.ileur  ealculee  pour  la  formation  a  1'equilibrc.  Lc  rapport  pour  Ic  carbone  cst  _ 0,62  ±  0  12'..  "C- 

jcllc  Va!j'"r  nCfii'"vc  c'*  aiiribnOo  a  la  formation  dc  methane  par  Jes  nombrcuscs  bacteries,  dc  telle  sorte  que  le’  s'ysteme 
i1  dont  0"  "a't  'I'"1  un  coefficient  dc  temperature  negatif,  dominc.  Lc  rapport  pour  Ihydrogenc  cst 

J.i  _  C,  environ  double  dc  cclui  que  Ton  peut  caleulcr  pour  la  formation  u  l'cquilibre  de  la  cellulose.  (Lcs  erreurs 

represement  une  deviation  standard). 
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ABSTRACT 


The  temperature  coefficient  of  three  stable  isotopes  in  a  blue-green  algae, 
growing  in  alkaline  hot  springs  at  temperatures  of  30°C  -  68°C,  have  been  1 

measured.  That  for  oxygen  is  found  to  be  +(2.2  ±1.0)  ppt/°C  which  is 
larger  than  the  value  computed  for  equilibrium  formation  as  is  to  be 
expected.  That  for  carbon  is  found  to  be  (-0.62  ±  0.12)ppt/°C;  the 
negative  value  is  attributed  to  methane  formation  by  copious  bacteria 
so  that  the  system  CH^  —  HCO^,  which  is  known  to  have  a  negative  tempera¬ 
ture  coefficient,  dominates.  That  for  hydrogen  is  found  to  be  (3.3  ±  0.7)ppt/°C, 
about  twice  that  computed  for  equilibrium  formation  of  cellulose.  (The  quoted 
errors  represent  one  standard  deviation.) 


INTRODUCTION 


Temperature  coefficients  of  some  stable  isotope  ratios  in  bio-organic 
material  have  been  measured  in  marine  plankton1’2*  in  peat3,  and  in 
wood  ,  and  are  listed  in  Table  I,  By  measuring  these  ratios  in  residual 
bio-organic  material  such  as  sediments,  peat  beds  and  tree  rings,  and 
combining  these  data  with  the  measured  temperature  coefficients  it  may 
be  possible  to  evaluate  temperatures  of  past  climates,  especially  if  the 
records  of  two  or  more  such  isotope  thermometers  are  evaluated  in  each 
data  base.  If  two  or  more  thermometers  record  like  changes  with  passage 
of  time  the  likelihood  is  enhanced  that  the  changes  are  truly  related  to 
changes  in  temperature.  A  second  point  is  that  the  temperature  coefficients 
which  have  been  measured  are  only  slightly  larger  than  those  computed3 
assuming  the  organic  material  to  have  been  manufactured  at  equilibrium 
(see  Table  II).  In  general,  the  observed  effects  may  be  expected  to  be 
larger  than  the  equilibrium  effect,  the  latter  being  the  ratio  of  the 
rates  for  the  forward  and  back  reactions.  In  reactions  occurring  far 
from  equilibrium,  only  the  forward  rate  applies. 

In  order  for  a  system  to  be  at  equilibrium  every  chemical  or  mechanical 
process  must  proceed  both  in  the  forward  and  reverse  directions  a  great 
many  times.  For  example,  if  a  reaction  proceeds  with  evolution  of  a  gas, 
and  if  the  gas  molecules  immediately  leave  the  site  of  formation  so  that 
the  revet se  reaction  cannot  occur,  there  can  never  be  equilibrium,  nor 
anv  reasonable  approximation  of  it.  Fractionation  of  isotopes  at  equili¬ 
brium  is  caused  by  the  difference  between  fractionation  by  the  forward 
rate  and  fractionation  by  the  reverse  rate,  the  two  rates  being  equal. 
Distillation  of  a  liquid  in  a  closed  container  is  an  example  of  fraction¬ 
ation  at  equilibrium.  Distillation  into  a  vacuum  is  kinetic  because  the 
reverse  reaction  is  completely  inhibited. 


The  temperature  coefficient  a  ,  of  an  isotope  ratio  R  is  defined  as  the 
derivative  of  R  with  respect  to  a  change  in  temperature  T,  AR/A  T  5  a. 


In  the  present  paper  we  report  measurements  of  the  temperature  coefficients 

1  fi  1  £ 

of  0  / 0  and  D/H  in  algae  growing  in  hot  springs  and  streams  at  Rotorua, 

4 

New  Zealand,  and  find  them  not  much  larger  than  that  measured  in  tree  wood  , 
(see  Table  1) . 

13  12 

The  ratio  C  /C  has  also  been  measured  but  its  temperature  coefficient 

4 

is  far  smaller  than  that  found  in  wood  ,  probably  because  the  decay  of 

13 

dead  algae  generates  methane  which  carries  off  C  leaving  depleted 
bicarbonate  to  feed  the  living  algae.  This  temperature  coefficient,  namely 
for  the  system  methane-bicarbonate,  reduces  the  overall  coefficient  to 
-0.62  ±  0.12  ppt/°C,  as  will  be  discussed  further  in  the  second  section 
of  this  paper. 

The  measurements  described  here  were  derived  from  a  limited  number  of 
samples  very  kindly  collected  by  Dr.  E.  F.  Lloyd^.  We  apologize  for  the 
small  number  of  the  data,  but  find  the  results  to  be  important,  even  if 
not  definitive,  because  they  l ear  out  the  fact  that  temperature  coefficients 
of  stable  isotopes  in  organic  material  are  measurable.  The  field  of 
isotope  thermometers  in  bio-organic  material  is  very  new,  so  that  each 
measurement  contributes  significantly  as  is  evident  in  Table  I. 

The  hot  springs  at  Rotorua  are  a  fascinating  bio-organic  system  which  will 
no  doubt  prove  rewarding  of  further  study. 
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EXPERIMENTAL 


Algae,  collected  at  several  places  in  the  hot  springs  at  Rotorua8,  New 

Zealand,  of  a  blue-green  variety^  were  freeze  dried  and  analyzed  mass 

18  16  13  12 

spectrometrically  for  0  / 0  and  C  /C  .  The  chemical  reactions 

Involved  are  as  follows: 

a.  For  measurement  of  0  /0  8,  HgCl2  and  algae  evolves  CO2  gas8  when 

heated;  HC1  is  removed  with  a  quinoline  trap. 

13  12 

b.  For  measurement  of  C  /C  ,  the  algae  were  burned  to  completion  in 
oxygen  gas,  yielding  C02. 

c.  For  measurement  of  D/H,  the  algae  were  heated  with  powdered  uranium 
metal  to  evolve  hydrogen  gas. 

For  each  measurement,  the  amount  of  reagent  used  was  much  larger  than 
needed  for  stoichiometric  proportions  with  5  mg  of  algae,  to 
insure  that  the  reaction  was  complete  so  that  there  would  be  no 
fractionation. 


The  measured  ratios  are  listed  in  Table  III  together  with  the  temperature 
of  the  water  in  which  the  algae  were  growing  when  collected.  The  isotope 
ratios  are  expressed  in  the  form, 


13 


(C13/C12) sample 

(c13/c12) 


-lx  1000 


’18 


(038/038) sample 

(018/016) 


-lx  1000 


(D/H) 


sample 


(D/H), 


-1  x  1000 


(1) 
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where  (0  /0  )Q,  (C  /C  fc)o  and  (D/H)q  have  been  chosen  as  the  values  of  the 

isotope  ratios  predicted  by  the  least  squares  ar.alysis^  for  50°C.  In  the 


analysis,  the  ratios  are  expressed  as  linear  functions  of  temperature, 
T,  according  to 


6  -  a  +  bT  (2) 

where  b  is  the  temperature  coefficient. 

The  values  of  a  and  b  computed  from  the  measurements  are  shown  in  Table  III. 

Since  the  water  reservoir  is  much  larger  than  the  CC^  reservoir,  it  is 
expected  that  most  of  the  oxygen  in  algae  will  be  supplied  from  the 
water.  It  is  known  for  Euglena,  a  water  plant,  that  this  is  the  case10. 

For  exchange  of  018  in  cellulose  with  018  in  water,  at  equilibrium,  the 
predicted  temperature  coefficient  is  +0.92  ppt/°C,  compared  to  which  the 
measured  temperature  coefficient,  +(2.2  ±1.0)  ppt/°C  reported  in  Table  I, 
is  positive  and  larger,  as  is  to  be  expected  for  a  natural  process  which 
is  likely  to  involve  dynamic  effects.  Since  oxygen  gas  is  evolved  in 
photosynthesis,  according  to  the  schematic  reaction, 

-°2  +  HjO  •  cellulose  +  0^  f 

and  leaves  the  site  of  the  reaction,  the  reverse  reaction  is  inhibited, 
and  equilibrium  considerations  would  not  be  applicable. 

The  temperature  coefficient  of  C13/C12  evaluated  from  the  measurements 
is  negative.  This  may  be  accounted  for  as  follows.  The  water  issuing 
from  the  hot  springs  is  strongly  alkaline  (pH~8.4)  and  in  addition 
contains  masses  of  decaying  algae  and  bacteria7  in  large  amounts,  so  that 
methane  is  being  evolved  by  decay  according  to  the  schema, 

HC03  +  bacteria  +  dead  cellulose  — ►  CH4  |  +  decaying  organic  matter 

Namely  the  system  is  dominated  by  exchange  of  C13  in  the  abundant  inorganic 
system 
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,13n_  ,  ,„12, 


HC~Or  +  (C^H. )  - ►  HC120-  r13 H  * 

3  4'dis«olved  «C  03  +  C  tydlMolvtd 


.13 


(C  ^dissolved 


(C13H.) 


V gaseous 


I 


so  that  the  liquid  system  is  depleted  of  C13. 


The  equilibrium  fractionation  factors  for  the  system  CC^  -  HCO-  -  CH  , 
calculated  by  Bottinga11,  show  that  C13  is  enriched  in  CH4  with  respect 
to  HC0-,  increasingly  at  high  temperatures,  and  that  the  minimum  value 
of  the  temperature  coefficient  (predicted  for  equilibrium)  is  negative, 
namely  -0.32  ppt/°C.  The  actual  temperature  coefficient  is  expected 
therefore  to  be  negative  and  somewhat  larger.  (Since  methane  is  evolving, 
dynamic  effects  clearly  are  playing  a  role.)  Thus  methane  evoluation  by 
bacterial  decomposition  of  the  dead  algae  carries  off  C13,  leaving  the 
aqueous  HC0~  depleted  in  C13  so  the  growing  algae  have  less  C13  at  high 
temperatures.  Farther  down  stream  where  the  water  is  cooler,  this 
process  of  enrichment,  being  temperature  dependent,  has  a  smaller  fraction¬ 
ation  factor.  Thus  in  this  case  the  temperature  dependence  of  C13/C12 

in  algae  acts  as  a  poor  thermometer,  one  appropriate  only  to  the  particular 
hot  springs  of  Rotorua,  New  Zealand. 
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PARTS  PER  THOUSAND 


Figure  1:  6^>  6Q  and  6^8  (relative  to  deltas  at  50°C)  vs.  water 

temperature. 


TABLE  1:  Experimental  values  of  temperature 

ratios  in  bio-organic  material. 


A.  C 


13 


B.  0 


Atmospheric  CC^ 

Wood  1.8  +  0.3 

Marine 

Phytoplankton  0.35 

Plankton  0.25 

Plankton  0.5 

Hot  Springs 

Bluegreen  Algae  -.62  +  .12 

18 


Atmospheric  Moisture 
Wood  A. A  +  0.3 

Hot  Springs 

Bluegreen  Algae  2.2  +  1.0 


C.  D 


Atmospheric  Moisture 
Wood  95  +  19 

Ground  Water* 

Peat  3 


Hot  Springs 


Blue  Green  Algae 


3.3  -  0.7 


*See  Reference  12  for  discussion  of  ground  water. 


coefficients  of  isotope 

Libby  and  Pandolfi 

Degens  et  al 
Sackett  et  al 

Eadie 

Libby  and  Pandolfi 

Libby  and  Pandolfi 

Libby  and  Pandolfi 

Libby  and  Pandolfi 

Schiegl 

Libby  and  Pandolfi 
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TABLE  II:  Temperature  coefficients  of  isotope  ratios  in  bio-organic 
material  computed  for  manufacture  at  equilibrium. 


COMPUTED*  FOR  CELLULOSE  ppt  1° C 
18 


a  <0  > 

0.96 

18 

a  (C02) 

1.14 

18 

a  (H20) 

0.92 

13 

a 

0.36 

a(HCOD) 

0.4 

a(OCOH) 

2.0 

♦FROM  REFERENCE  5. 
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TABLE  III:  Delta  values  for  isotope  ratios  in  blue  green  algae  growing 
in  hot  springs,  Rotorua,  New  Zealand,  and  temperature 
coefficients  evaluated  from  the  measurements. 


T 

degrees 

Centigrade 

S13* 

parts  per 
thousand 

S18* 

parts  per 
thousand 

D 

parts  per 
thousand 

30° 

+13.7 

— 

-56 

46 

-2.6 

-23.3 

-29  ! 

51 

+1.4 

+17.0 

-13 

55 

-1.2 

+15.1 

-9 

56 

-3.3 

+17.7 

+22 

64 

-9.3 

+22.4 

+64 

*  The  deltas  are  the  differences  from  the  values  predicted 

o  9 

for  50  C  by  least  squares  analysis  . 


S13  ' 

[31.1  -  [0.62  t  0.12] 

S18 

[  -110  +  [2.2  !  1.0]  T 

;  PPt 

6d  s 

[-165  +  (3.30  t  0 . 7 3 ) t} 

ppt 

The  relatively  large  errors  in  the  evaluated 

ppt 


relate  for  the  most  part  to  fluctuations  in  water  temperature  issuing  from 
the  hot  springs  as  described  in  Reference  7. 
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abstract 


Sr*, te  t  “”"ai  c"m“' 

records  for  these  years,  from  England  nearby  Basel  and  Henevi  i  ,h-54  AlD’  ,3nd  correIa,ei1  w»h  «hc  existing  wcathei 
arc  evaluated  as  (5.29  *  0.68),  (2.91  *  Odl/and  '  1  th.s  way>  the  temperature  coefficients  for  0'*/Qi* 

as  (2.73  +  0.67),  (2.01  ^  0.37)  and  “37  + n  am  l  *  d'52)  r«PecllVeIy.  the  coefficients  for  OVC»  arc  evaluated 

units  of  parts  per  thousand  per  degree*  Centigrade' ' (pptrCj'.'Vse  o7  these  th^’  **  189  It 16)*  •  167  ~  6)  and  171  ±  8)  in 
may  allow  the  study  of  air  temperature  changes  in  climate  before  1712.  *  m‘mcncaI,y  caI‘brated  isotope  thermometers 


R£sum£ 


at  Chuquc  tSJZS  de  V^sTa'd  T?  ?“ 

**  *  ^  vary: 

'««■“*»•*  P0“r  mC/.-'C  sont'respcaivemen,  ^??9  ^?  ?0?PoT\,2?xi°*Zi  2‘91  *0-41  «  2.86  ±  0.52 

(•7  ±  6  et  71-8  en  pour  rnille  par  "C.  L’utilisation  ’  de  ces  thermonZir'*7  ■  r37  *  <*41  et  pour  D/H89’±J6 

permettre  letude  des  variation,  de  la  temperature  de  Pair  avant  1712  tsotop.ques  numeriquement  etalonnes  peui 


tcntrale. 
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tn  a  journal  of  very  broad  scope,  to  avoid  the  use  of  specialized  jargon  that  would  be  understood  only  bv  a  small  group  of  worke-s 
in  a  narrow  field.) .  '  6  r 
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ibUlUft  THERMOMETERS  IN  OAK  TREE  RINGS  CALIBRATED 
USING  OFFICIAL  WEATHER  RECORDS 


Annual  tree  rings  are  a  most  valuable  source  of  information,  perhaos 
tne  only  year-to-year  source  of  information  on  environmental  and  climatic  changes 
during  tne  past.  There  are  now  tree-ring  sequences  available  that  go  back  more' 
than  7,000  years,  and  the  study  of  this  material  is  most  desirable.  ~v=-y 
effort  to  use  it  is  most  commendable.  This  paper  reports  the  d'seove-y  or 
an  impressive  correlation  between  stable  isotope  ratios  in  the  wood  of  an  oak 
and  temperature  records  in  Europe  for  the  18th  and  19th  centuries.  Should  this 
correlation  turn  out  to  be  a  general  phenomenon  observable  in  any  kind  of  wood, 
tren  it  would  provide  the  best  possible  way  for  determining  climatic  variations 
ror  the  time  prior  to  that  when  meteorological  records  started.  I  strongly  rec¬ 
ommend  chat  this  paper  be  published,  as  it  undoubtedly  will  sf‘-nil=>t=*  -urth°r 
research^n  this  important  field,  which  is  of  interest  in  a  lerv9  number  o’f~ 

fields  oi  science,  i  consider  the  length  of  the  paper  and  the  wav  the  findings 
are  presented  to  oe  adequate.  8 
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ABSTRACT 


The  stable  isotope  ratios  of  carbon,  oxygen  and  hydrogen  have  been  measured 
for  an  oak  from  central  Germany  in  wood  samples  consisting  of  roughly  three 
years  each,  for  the  years  1712  -  1954  A.D.  and  correlated  with  the  exist¬ 
ing  weather  records  for  these  years,  from  England,  nearby  Basel,  and 
Geneva.  In  this  way,  the  empirical  temperature  correlations  can  be 

1  O  If. 

expressed  by  the  following  coefficients  for  0/0:  (5.29  +  0.68), 

(2.91  +  0.41)  and  (2.86  +  0.52)  respectively,  the  corresponding  coeffi¬ 
cients  for  C13/C12  are  (2.73  +  0.67),  (2.01  +  0.37)  and  (2.37  +  0.41), 
respectively,  and  for  D/ll  (89  +  16),  (67  +  6),  and  (71  +  8)  in  units  of 
parts  per  thousand  per  degree  Centigrade  (ppt/°C).  Use  of  these  three 
numerically  calibrated  Isotope  thermometers  may  allow  the  study  of  air 
temperature  changes  in  climate  before  1712. 


INTRODUCTION 


Information  on  climate  is  stored  in  natural  data  banks,  namely  in  bio-organic 
materials  that  are  deposited  year  after  year  and  thereafter  remain  undisturbed. 
In  particular,  tree  rings  in  living  and  dead  trees  may  provide  a  worldwide 
record  of  climate  for  the  last  few  thousand  years.  The  information  is  in  chemi¬ 
cal  form,  namely  in  the  ratios  of  stable  isotopes;  that  is,  variations  in  air 
temperature,  and  therefore  in  the  temperature  at  which  the  wood  is  formed, 
cause  variations  in  the  stable  isotope  ratios  in  the  new  wood  as  it  grows. 

An  additional  variable  is  the  temperature  of  formation  of  water  which  nourishes 
the  tree. 

The  justification  for  regarding  stable  isotope  ratios  in  land  plants  as  thermom¬ 
eters  for  air  temperature  is  twofold.  Although  we  also  need  to  know  the  isotope 
ratios  in  atmospheric  C02  and  of  rainwater  which  nourishes  the  plants, there  is 
good  reason  to  believe  that  in  the  long  term  average,  atmospheric  CC>2  has  re¬ 
mained  unchanged  for  a  long  time,  say  millions  of  years,  so  this  is  not  an 
unknown.  The  isotope  ratios  in  precipitation  have  been  studied  all  over  the 
world  s  surface  for  the  last  15  years  and  have  been  reported  in  the  data  compila¬ 
tions  of  the  International  Atomic  Energy  Agency.  So  we  know  the  percent  of 
18  16 

0  / 0  and  D/H  in  rainwater  as  a  f-n«.  tion  of  altitude  and  latitude  and  longi¬ 

tude  of  the  surface  of  the  earth.  These  ratios  may  have  undergone  variations  in 
the  past  caused  by  changes  in  the  global  temperature.^" 

Consequently,  stable  isotope  ratios  in  land  plants  depend  on  several  temperature 

dependent  functions.  The  resulting  temperature  dependence  has  been  observed 
experimentally. 

When  a  plant  dies,  the  record  of  the  air  temperature  may  be  stored  in  its 
tissues,  in  the  form  of  the  numerical  values  of  the  stable  isotope  ratios. 
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The  most  important  of  these  ratios  are  D/H,  C13/C12,  018/016,  N15/N14,  and 
S34/332. 

The  ratios  were  measured  for  three  independent  elements  in  the  wood,  H,  C,  0, 
that  is,  the  records  from  three  thermometers  were  read,  so  that  it  appears 
that  temperature  variations  in  the  past  can  be  evaluated,  with  some  confidence 
and  the  observed  variations  seem  truly  to  be  correlated  with  temperature  changes. 

Factors  that  can  be  calculated  quantitatively  are  the  equilibrium  temperature 
coefficients  of  the  stable  isotopes  of  H,  C,  and  0  in  cellulose  (trees  being  69% 
cellulose,  the  remainder  being  mainly  lignin).  They  are  large  enougi.  so  that 
variations  of  a  few  degrees  in  air  temperature  should  cause  measurable  effects.2 
The  values  of  the  coefficients  for  D/H,  C13/C12,  and  018/016  have  been  computed 
in  Reference  2  (see  Table  1)  on  the  assumption  that  cellulose  is  formed  at 
chemical  equilibrium.  It  is  known  that  at  least  for  C13/C12  in  marine  plankton3*4, 
and  for  D/H  m  peat,5  the  measured  coefficients  are  approximately  equal  to  the 
computed  values.  In  general,  however,  the  observed  effects  may  be  expected  to 
be  larger  than  the  equilibrium  effect,  the  latter  depending  on  the  difference 
of  the  rates  for  the  forward  and  back  reactions.  Two  temperature  effects  are 
probably  involved,  the  first  characterizing  manufacture  of  cellulose  and  the 
second  characterizing  manufacture  of  rain  water.  Eventually  these  effects 

should  be  separated.  Until  now,  the  temperature  coefficients  in  trees  had  not 
yet  been  measured. 

Anyhow,  because  of  variations  in  018/016  and  D/H  rati  in  precipitation,  the 
stable  isotope  thermometers  in  trees  must  be  calibrated;  namely,  their  observed 
variations  in  rings  formed  in  the  last  few  hundred  years  must  be  compared  with  the 
official  weather  records  of  variations  in  air  temoerature  in  !he  neighborhood  of 
the  tree,  to  yield  numerical  values  for  the  changes  in  isotope  ratios  per  degree 
Centigrade  change  in  temperature. 
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Calibration  of  the  effects  of  climate  on  stable  isotope  ratios  in  living  trees 
must  be  undertaken  judiciously,  keeping  in  mind  that  numerical  air  temperatures 
have  been  recorded  at  best  for  no  more  than  about  three  hundred  years  in  Europe, 
and  for  only  about  one  hundred  years  in  America.  Furthermore,  both  the  tempera¬ 
ture  records  and  the  recently  formed  tree  rings  may  be  perturbed  by  local  climate 
changes  caused  by  increasing  heat  and  air  polli  tion  produced  by  the  cities  growing 

up  around  the  existing  weather  stations  and  around  the  trees.  Also  the  introduction 

13  12 

of  fossil  fuel  CO2  (Suess  Effect)  is  expected  to  perturb  the  C  /C  ratios  in 
rings  grown  after  1890. 

13  12  18  16 

In  the  present  paper  are  reported  measurements  of  C  /C  ,  0  /0  ,  and  D/H  in 

the  rings  of  an  oak^  which  grew  in  Aalen,  Germany,  10°11'  east,  48°50'  north, 
about  half  way  between  Mannheim  and  Munich,  from  1700  to  1965  A.D.  Air  tempera¬ 
ture  records  exist  for  nearby  Basel ^  since  1755,  for  Ceneva^  since  1768,  and  for 

g 

England  from  1698.  The  measured  isotope  ratios  as  a  function  of  tree  ring  data 
have  been  compared  with  these  official  temperature  records  for  1712  -  1954  A.D. 
using  the  dated  tree  ring  sequence^  to  determine  the  corresponding  age  of  the  tree 
rings  in  years.  In  this  way,  temperature  coefficients  of  the  stable  isotope 
ratios  have  been  evaluated,  as  a  phenomenological  effect,  without  attempting 
to  separate  the  effect  of  manufacture  of  cellulose  from  that  for  manufacture 


of  precipitation. 


THEORETICAL 


Of  the  various  falters  that  determine  the  isotope  ratio  for  a  given  element 
in  plant  material,  at  least  one,  the  thermodynamic  isotope  exchange  constant, 
can  be  calculated  rigorously.  Consider,  for  example,  exchange  of  C02  and 
cellulose  at  equilibrium  according  to  the  schematic  reaction: 

Consider,  for  example,  exchange  of  C02  and  cellulose  at  equilibrium  according 
to  the  schematic  reaction: 

C1302  +  (H-C12  -  0H)n - *C1202  +  <HC13  -  0H)n 

The  equilibrium  constant  K  is  given  by, 

K  -  _Q(HC130H) 

Q(c13o2)  q(hc12oh) 

where  the  0's  depend  only  on  temperature  according  to 

-U/2)tWkT 

Q  •=  _ _ 

[1  .  e-Wklj 

for  each  frequency. 

These  w’s  are  frequencies  of  vibration,  and  depend  on  the  force  constant  k  and 
the  reduced  mass,  m,  of  each  two  vibrating  atoms,  according  to: 

U)  -  2 

m 

For  a  given  pair  of  atoms  bound  to  each  other,  the  force  constant,  or  spring 
constant  k,  has  the  same  value  regardless  of  isctopic  mass.  The  isotopic  effect 
occurs  because  m  is  different  for  different  masses,  so  the  m's  are  different. 
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For  the  C-0  bond  for  example,  there  are  four  different  frequencies  corresponding 
to  the  isotopic  pairs, 


12n16  c13016  „12n18  „13„18 


C  O 


C  O 


C  O 


Thus,  in  land  plants  the  equilibrium  constant  K  depends  on  air  temperature,  T. 
The  equilibrium  constants  K(300°C)  and  K(273°C)  for  the  formation  of  cellulose 
at  equilibrium  have  been  computed  according  to  the  schematic  reaction: 

C02  +  H20  —  (H  -  C  -  0H)n  +  02 


The  temperature  coefficient  a  for  exchanges  of  the  stable  isotopes  between 
cellulose  and  C02>  or  H20,  or  02  is  then  obtained  from  the  computed  K's  according 
to: 

a  -  K(300°C)  -  K(273°C) 

27°C 

2 

These  computed  coefficients  for  exchange  at  equilibrium  are  found  to  be  large 
enough  to  measure,  see  Table  1.  If  exchange  occurs  away  from  equilibrium,  (these 
are  called  kinetic  effects),  the  temperature  coefficients  can  be  larger. 

An  example  of  a  kinetic  effect  is  evolution  of  a  gas,  e.g.,  0 2  in  the  above 
example,  preventing  the  reverse  reaction  from  occurring  and  hence  preventing 
equilibrium  conditions  to  be  reached. 

In  an  equilibrium  process,  atoms  and  molecules  must  react  and  back-react  many 
times.  If  a  substance  leaves  the  site  of  the  reaction  as  soon  as  it  is  formed, 
there  can  never  be  equilibrium.  For  instance,  if  water  vapor  escapes  into  a  vacuum 
at  the  instant  it  leaves  the  surface  of  liquid  water,  or  if  oxygen  leaves  a  tree 

leaf  as  soon  as  formed  in  photosynthesis,  these  are  dynamic  effects  and  preclude 
equilibrium. 
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Of  the  five  principle  isotope  thermometers  in  bio-organic  material,  no  two 
atoms,  H,  C,  0,  N  and  S  are  tied  through  the  same  set  of  bonds.  For  example, 
in  cellulose  H  is  tied  to  C  and  to  0;  0  is  tied  to  H  and  C;  and  C  is  tied  to  0 
and  H  and  other  carbons.  Thus  there  are  three  independent  thermometers,  and 
two  more  if  N  and  S  are  added.  In  shell  there  is  only  one  thermometer  because 
C  and  0  are  tied  through  the  same  bond.  If  one  isotope  ratio  changes,  there 
is  uncertainty  as  to  the  cause,  hut  if  three  isotope  ratios  depend  on  temperature 
as  well  as  on  other  parameters,  and  all  three  change  simultaneously,  then  it  is 
likely  that  the  changes  were  caused  by  temperature  change. 
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EXPERIMENTAL 


In  a  2"  thick  slice  cut  across  the  tree,  a  groove  1/3"  wide  and  1/3"  deep  was 
milled  perpendicular  to  the  tree  rings,  that  is  along  a  radius  of  the  tree,  from 
the  innermost  (oldest)  part  of  the  tree  out  to  the  bark  (see  Figure  1) .  Each 
time  that  the  milling  tool  had  moved  1/3",  the  sawdust  produced  was  collected 
into  an  individual  vial  with  the  aid  of  a  camel’s  hair  brush.  Each  sample 
consisted  of  sawdust  from  about  3  or  4  rings,  depending  on  the  variations  in  ring 
widths,  the  total  number  of  samples  being  68  for  the  time  span  1712  -  1954  A.D. 

The  stable  isotope  ratios  were  measured  in  the  sawdust  samples  by  mass  spectro¬ 
metry.  The  spectrometer  has  a  very  small  quartz  manifold,  of  volume  about 
3 

100  cm  .  About  5  mg  of  sawdust  were  used  for  each  measurement.  The  chemical 
. eactions  used  were  as  follows: 

a.  For  measurement  of  C13/C12,  the  sawdust  was  burned  to  completion 
in  oxygen  gas,  yielding  C02 . 

b.  For  measurement  of  0  ^/0^6,  HgClj  +  sawdust  evolves  CO^  gas^  when 
heated;  HC1  is  removed  with  a  quinoline  trap. 

c.  For  measurement  of  D/H,  sawdust  is  heated  with  powdered  uranium, 
to  evolve  Hj. 

For  each  measurement,  the  amount  of  reagent  used  was  much  larger  than  needed  for 
5  mg  of  sawdust,  to  insure  that  the  reaction  was  complete  so  that  there  would  be 
no  fractionation.  The  residue  was  inspected,  after  cooling,  for  the  same  reason. 
The  measured  ratios  are  listed  in  Table  I  along  with  the  tree  ring  dates  as 
determined  by  the  Munich  Forest  Botanical  Institute6.  The  measured  variations 
in  stable  isotope  ratios  in  tree  rings  are  seen  to  be  large,  of  the  order  of  a 
few  parts  per  thousand  (ppt) ,  expressed  in  the  form: 
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x  1000 
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(C13/C12)  sample 
(C13/C12)q 


6 


18 


6 


D 


(0l8/018)  sample 

(018/016)q 

(D/H)  sample 
(D/H)q 


x  1000 


x  1000 


The  normalizing  ratios  (C13/C12)0,  (018/016)Q  and  (D/H)q  may  be. 
ient  way,  and  here  have  been  cnosen  arbitrarily  as  the  measured 
years  1712  -  1714  A.D.  for  computing  the  delta  values  listed  in 


defined  in  any  conven- 
ratios  for  the 
Table  I. 


The  snpwood  of  the  oak  extended  over  about  30  tree  rings.  That  is,  sapwood 
changed  into  inert  hardwood  so.i.e  30  years  after  it  was  manufactured,  and  thus 
for  about  30  years  each  tree  ring  was  bathed  in  flowing  sap  from  ground  water. 
Consequently,  it  is  to  be  expected  that  oxygen  in  the  alcohol  groups  of  cellulose, 
its  modular  formula  being  (HCOH^,  is  exchanging  with  sap  water,  and  thus  that 
the  stable  isotope  ratio  of  oxygen  in  each  ring  is  averaged  over  about  30  years1? 
There  is  less  reason  to  believe  that  carbon  in  the  wood  is  exchanging  with  carbon 
dioxide  dissolved  in  sap,  but  it  has  not  been  proven  that  there  is  no  exchange11,12. 
On  this  account,  the  measured  isotope  ratios  are  reported  in  the  form  of  a  running 
average  over  9  samples  at  a  time,  namely  as  an  average  over  the  time  that  a 
given  ring  is  sapwood.  Winter  air  termparatures,  from  official  weather  records7,8 
averaged  in  the  same  way,  are  also  listed  in  Table  1.  Temperatures  for  January, 
February  and  March  show  the  major  increase  since  the  Little  Ice  Age  (1440-1850  A.D.), 
whereas  the  summer  temperatures  then  were  not  much  different  from  now13.  So  the 
average  yearly  temperatures  vary  in  the  same  way  as  the  average  temperatures  of 
January,  February  and  March,  but  with  a  somewhat  smaller  amplitude.  Oak  trees  make 
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the  greatest  growth  of  the  year  in  early  spring  when  the  ground  thaws  and  the 
winter's  accumulation  of  precipitation  melts,  so  it  seems  reasonable  that  they 
will  be  sensitive  to  late  winter  temperatures,  especially  for  isotopes  of  oxygen 
and  hydrogen  which  come  from  the  melt  water  and  depend  on  the  temperature  at  which 
precipitation  formed. 

The  carbon  ratios  (6.^)  for  the  years  1890  -  1950  A.D.  have  been  corrected 

for  fossil  carbon  production,  the  maximum  correction,  that  for  the  year  1950  A.D., 

taken  as  +8.4%  or  a  maximum  increase  of  +2.1  ppt  in  6^*  This  correction  was 

estimated  as  follows.  In  wood  from  rings  of  1920  A.D.  two  radiocarbon  dates 
14  + 

were  measured  as  375  _  35  years  old  with  respect  to  1950  A.D.  ,  whereas  the 

actual  age  is  1950  -  1920  =  30  years,  corresponding  to  4.2%  dilution  of  atmospheric 
14 

C  0^  by  inert  CO^  produced  by  man's  burning  of  coal  and  oil  up  to  1920.  In 

1950  the  correction  for  at  this  particular  place  should  be  8.4%  and  the 
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correction  for  C  dilution  in  atmospheric  therefore  should  be  8.4%  of 
14 

25  ppt  or  2.1  ppt;  see  further  discussion  in  Reference  15. 

The  ratios  and  temperatures  were  correlated  by  least  squares  analysis,  expressing 
the  ratios  in  the  form, 

6  =  a  +  bT  [2] 

and  evaluating  the  coefficients  a  and  b,  (a  is  the  value  of  6  at  T  =  0,  and 

b  is  the  temperature  coefficient).  Their  numerical  values  are  shown  in  Table  II 

and  plotted  in  Figures  2  and  3.  The  average  of  the  temperature  coefficients 
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for  0  /0  computed  from  air  temperatures  at  Basel  and  Geneva,  (these  places 

being  closer  to  the  tree  than  is  England)  is  (2.88  t  0.36)ppt/°C.  This  may  be 
compared  to  the  value  of  2.2  t  1.0  ppt/°C  measured  for  blue-green  algae^.  The 
average  for  Ci7C  in  oak  is  2.19  -  0.25  ppt/°C,  which  is  a  good  deal  larger 
than  the  value  of  0.35  ppt/°C  measured  by  Eadie^  for  marine  phytoplankton. 
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The  oak  tree  oxygen  depends  on  changing  rain  water  however,  whereas  the  0  and  H 
in  the  water  plants  should  be  from  an  unchanging  water  supply.  The  average  for 
D/H  is  69  +  4,  and  it  also  is  large  compared  with  the  value  of  3  ppt/°C  derived 
from  Schiegl's  measurement.  A  comparison  is  made  of  the  averages  of  coefficients 
for  England,  Basel,  and  Geneva  temperatures  with  theoretical  coefficients  in 
Table  III.  In  a  manner  of  speaking,  trees  store  old  rainwater.  Co,  to  the  isotope 
fractionation  caused  by  evaporation  of  water  vapor  from  the  oceans  and  precipita¬ 
tion  as  rain,  the  tree  adds  a  further  fractionation  caused  by  absorption  of  ground 
water  into  the  salt  and  sugar  solutions  in  the  root  system,  and  by  manufacture  into 
bio-organic  material.  The  combination  of  these  processes  appears  to  produce  large 
temperature  coefficients. 

If  the  correlations  reported  for  this  oak  are  general  for  all  oaks,  then  this 
effort  can  be  very  useful. 
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YEARLY  FLUCTUATIONS,  STftNAT.  "NHTQF" 


Yearly  variations  of  C13/C12  in  atmospheric  C02  resemble  a  sine  function  with  a 
few  parts  per  thousand  peak  to  peak  amplitude;  the  amplitude  may  be  geographically 
dependent17.  All  available  data  on  this  seasonal  variation  appear  to  show 
depletion  of  C  in  the  winter,  enhancement  in  summer,  and  the  greatest  rate  of 
change  in  spring  and  fall.  Thus,  if  spring  comes  early  or  late,  namely  if  the 
tree  begins  it  major  growth  of  a  new  ring  early  or  late  with  respect  to  the  inflection 

point  of  the  sine  curve,  the  C13/C12  ratio  in  the  new  wood  may  be  depleted  or  enhance- 
cd  by  up  to  3  ppt,  with  respect  to  the  average  ratio. 

One  expects  that  similar  yearly  variations  occur  for  018/016  in  atmospheric  C02. 

For  comparison  it  is  known  that  there  is  a  similar  seasonal  variation18  both  for 
C02  /C0,K  and  for  C1302/C1202  in  the  atmosphere  of  a  few  parte  per  thousand. 

However,  this  source  of  year  to  year  variations  in  the  oxygen  isotope  ratio  in 
wood  is  probably  less  than  the  variations  in  rain  and  snow. 

Namely,  at  any  given  geographical  place,  the  018/016  ratio  in  rain  water  varies  by 
some  10  ppt  during  the  yea^9  ,  even  fluctuating  by  several  parts  per  thousand  month 
by  month  and  day  by  day,  depending  on  the  temperatures  at  which  sea  water  evaporated 
and  at  which  the  precipitation  formed,  and  therefore  on  the  altitude  of  the  clouds 
(see  Figure  4).  This  dependence  on  temperature  is  exact  but  the  dependence  on  air 
temperature  at  ground  level  is  valid  only  on  the  average  (see  Figure  5),  indicating 
that  the  average  air  temperatures  aloft  where  precipitation  forms  are  correlated  to 
average  air  temperatures  at  ground  level.  At  higher  latitudes  where  trees  make 
a  rapid  spring  growth  consuming  the  melt  water  from  accumulated  winter  snow, 
both  that  o-  the  surface  of  the  ground  and  that  stored  in  the  soil  as  winter  frost, 
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the  new  tree  ring  binds  oxygen  whose  isotope  ratio  represents  an  average  over 
the  winter  time  precipitation. 

If  the  winter  has  been  colder  than  average,  the  melt  water  is  depleted  in  018, 

and  if  warmer  than  average,  then  it  is  enriched,  to  the  degree  of  several  parts 

per  thousand  (see  Figure  4  and  5).  Study  of  the  IAEA  data  shows  that  5  is 

18 

independent  of  total  rainfall. 

It  is  because  of  the  large  yearly  fluctuations  in  the  air  temperature  data  that  we 
have  plotted  their  running  averages  and  used  them  in  the  least  squares  fits. 

Comparison  of  a  plot  of  8lg  vs  in  tree  rings  (see  Figure  6)  with  that  in 
precipitation  'hows  that  the  slope  of  8  for  precipitation  no  longer  obtains, 
probably  because  the  rates  of  exchange  of  D  and  018  between  cellulose  and  sap 
are  quite  different  from  each  other.  (Thus,  it  may  be  now  possible  to  study 
dynamic  physiology  of  trees,  as  a  research  independent  of  the  study  of  climate.) 

In  the  long  term  of  tens  of  hundreds  of  years,  a  distinct  climate  change  may 
record  itself  by  a  distinct  change  in  the  018/016  ratio  in  precipitation  and 
therefore  in  the  tree  wood  manufactured  after  the  climate  change  compared  with 
that  from  before.  It  is  known  that  in  a  water  plant,  Euglena,  90%  of  the  bound 
oxygen  is  suppl  from  water  .  In  trees,  the  source  of  oxygen  is  not  known, 
but  because  the  amount  of  water  in  the  tree  is  so  much  greater  than  the  amount 
of  dissolved  carbon  dioxide,  it  is  to  be  expected  that  the  major  source  is 
rainwater,  even  if  CO2  and  I^O  are  in  equilibrium. 
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CONCLUSION 


Every  kind  of  thermometer  has  to  be  calibrated.  The  stable  isotope  ratios  in 
tree  wood  depend  partly  directly,  partly  indirectly  on  the  temperature  at  which 
the  wood  was  manufactured.  The  present  paper  reports  the  phenomenological  cali¬ 
bration  of  the  oxygen,  carbon  and  hydrogen  isotope  ratios  in  a  European  oak.  The 
long  term  increase  in  air  temperatures  from  the  Little  Ice  Age  to  the  present 
offered  the  opportunity  to  make  this  calibration.  In  each  case  the  measured 
temperature  coefficients  are  larger  than  those  computed  assuming  the  wood  to  have 
been  manufactured  in  equilibrium  with  C02  and  H^,  as  is  to  be  expected  because 
kinetic  effects,  no  doubt,  play  a  role,  and  because  in  formation  of  precipitation 
used  in  the  manufacture  of  cellulose  the  temperature  coefficients  of  018  and  D 
are  large  and  positive.  In  this  first  attempt  to  study  isotope  thermometers  in 
tree  rings  we  have  lumped  the  temperature  coefficients  for  cellulose  manufacture 
with  the  temperature  coefficients  for  precipitate  formation  for  an  evaluation  of 
the  over-all  phenomenological  coefficient.  We  realize  that  many  problems  must 
be  solved  for  a  better  understanding  of  these  effects,  that  woody  material  must 
be  grown  under  controlled  conditions  for  evaluation  of  the  cellulose  manufacture 
coefficients,  that  other  oaks  must  be  compared  with  this  oak,  and  that  oaks 
should  be  comparec  with  other  kinds  of  trees.  Using  the  measured  temperature 
coefficients  of  these  thermometers,  and  using  measured  isotope  ratios  in  chrono¬ 
logic  sequences  of  tree  rings  manufactured  before  1712,  it  now  may  be  possible  to 

evaluate  the  European  climate  from  Lhe  centuries  before  the  advent  of  mechanical 
22 

thermometers  .  Use  of  more  than  one  thermometer  increases  the  likelihood  that 
temperature  effects  are  being  observed.  Other  isotope  thermometers  in  tree  rings, 
namely  nitrogen  and  sulfur,  should  also  have  raeasureable  temperature  coefficients. 
The  authors  are  deeply  grateful  for  helpful  suggestions  from  Samuel  Epstein, 

Hans  Suess,  Ian  Kaplan,  Harmon  Craig,  Devendra  Lai,  Edward  Teller  and  Willard  Libby. 
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TABLE  I 


Measured  isotope  ratios  in  dated  oak  tree  rings  and 
corresponding  winter  temperatures  from  official  weather 
records  for  England,  Basel  and  Geneva. 


SAMPLE 

NUMBER 

613 

ppt 

60 

ppt 

618 

ppt 

TEng°C 

T 

'Basel  L 

^Geneva 

cC 

DATE  SPAN 

1 

0 

0 

0 

4.26 

1712  -  1714 

2 

1.13 

156.0 

3.62 

3.89 

1715  -  1717 

3 

1.90 

181.0 

6.35 

3.98 

1718  -  1720 

4 

3.41 

141.0 

6.35 

4.24 

1721  -  1723 

5 

3.56 

110.0 

7.20 

4.36 

1724  -  1726 

6 

3.87 

131.0 

7.96 

4.30 

1727  -  1731 

7 

4.32 

103.0 

8.36 

4.39 

1732  -  1734 

8 

4.02 

73.2 

8.28 

4.46 

1735  -  1737 

9 

4.09 

51.1 

7.88 

4.38 

1738  -  1743 

10 

4.55 

23.2 

7.12 

4.52 

1744  -  1747 

11 

4.47 

46.4 

6.67 

4.52 

1748  -  1750 

12 

i  3.49 

31.7 

7.24 

4.25 

1751  -  1753 

13 

3.71 

58.0 

5.39 

4.00 

0.06 

(0.14) 

1754  -  1758 

14 

3.11 

63.5 

5.59 

4.03 

0.98 

(2.20) 

1759  -  1761 

15 

3.03 

45.7 

5.11 

3.98 

0.81 

(1.90) 

1762  -  1764 

16 

3.41 

54.2 

4.38 

3.95 

0.93 

(2.30) 

1765  -  1767 

17 

4.32 

60.3 

4.18 

3.97 

1.00 

2.34 

1768  -  1772 

18 

55.7 

3.29 

3.86 

0.95 

2.72 

1773  -  1775 

19 

4.40 

62.7 

4.50 

3.73 

1.02 

2.49 

1776  -  1777 

20 

3.64 

68.9 

3.21 

3.70 

1.21 

2.25 

1778  -  1779 

21 

4.24 

75.1 

2.09 

3.69 

IB 

11 

1780  -  1783 

22 

4.85 

48.0 

2.01 

3.84 

H3< 

na 

1784  -  1786 

23 

5.46 

41.8 

1.93 

3.77 

1.12 

2.17 

1787  -  1791 

24 

5.46 

37.9 

1.85 

3.82 

1.20 

2.00 

1792  -  1794 

25 

3.49 

39.5 

3.78 

3.82 

1.33 

2.07 

1795  -  1756 

26 

3.00 

29.4 

3.94 

3.66 

HI 

1.85 

1797  -  1801 

27 

2.73 

51.1 

4.22 

3.86 

o 

2.09 

1802  -  1804 

28 

2.12 

45.7 

3.78 

3.89 

WBM 

2.11 

1805  -  1808 

29 

1.21 

15.5 

3.22 

3.93 

;TO 

1.90 

1809  -  1811 

30 

1.36 

10.1 

4.46 

3.99 

TO 

1.98 

1812  -  1814 

TABLE  I  (continued) 


TABLE  I  (continued) 


SAMPLE 

NUMBER 

H 

6D 

PPt 

618 

PPt 

T  °C 
'Eng  C 

T  °C 

'Basel  u 

ESMj 

date  span 

61 

7.05 

B3S 

9.69 

mm 

2.29 

2.95 

62 

7.05 

9.37 

■9 

2.10 

2.85 

63 

7.05 

155.0 

9.12 

4.42 

2.10 

2.92 

1939  -  1941 

64 

7.20 

159.0 

9.32 

4.43 

2.17 

2.95 

1942  -  1943 

65 

7.50 

159.0 

9.33 

4.41 

2.38 

3.14 

1944  -  1945 

66 

7.66 

154.0 

7.20 

4.41 

2.50 

2.98 

1946  -  1947 

67 

7.96 

156.0 

6.15 

4.30 

2.55 

3.34 

1948  -  1950 

68 

10.31 

108.0 

9.89 

4.24 

2.19 

2.95 

1951  -  1954 

Average  standard  deviation  on  each  9  sample  average  of  lg  =  0.31  ppt 
Average  standard  deviation  on  each  9  sample  average  of  «  0.21  ppt 

I  0 


TABLE  I 

Nine  sample  running  averages  for  measured  and  ,Q  vs. 

year  of  center  of  the  sample  and  vs.  official  average  winter 
temperatures  (Jan.,  Feb.,  March)  in  England  and  Basel  and 
Geneva.  Temperatures  for  Geneva  for  1759  -  1783  (shown 
in  parentheses)  were  not  recorded  and  have  been  chosen  by 
normalizing  to  Basel  temperatures. 
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TABLE  II 


Temperature,  coefficients  computed  from  measured  isotope 
ratios  and  weather  records. 


EXPERIMENTAL  TEMPERATURE  COEFFICIENTS** 

Computed  by  correlating*  the  isotope  ratios  with  official6,7  average 
winter  temperatures  (averaged  over  January,  February  and  March). 

ISOTOPE 

RATIO 

TEMPERATURE  COEFFICIENT 

Number  of 

Samples 

68 

018 

o16 

*  t(5.29  t  0.68)TEng,and]ppt 

+  C(291  t  0.41)TBasel]ppt 

56 

♦  [(2.86  i  O.S2)TeerievaJppt 

56 

c13 

c12 

*  [(2-23  -  0.67)TE„g1and]ppt 

68 

+  C(2.01  +  0.37)TBasel]ppt 

56 

*  [  (2-37  i  0.41)TCeneva]ppt 

56 

D 

H 

*  [(89.5  t  16.0)TEng1andJppt 

68 

♦  [  (67.4  i  6.4)TBase,]ppt 

56 

*[  171  •‘l  ; 

56 

Temperature  Coefficients  in  (ppt)  •  (°C)'1 


*The  Bio  Medical  Data  Processing  Program20  obtained  from  the  UCB  computer 
department,  courtesy  of  Dr.  James  Enstrom,  was  used  to  correlate  the  isotope 
ratio  measurements  with  the  official  air  temperatures. 

**I^Ltemperatf-  c?efficie"ts  were  evaluated  using  9-sample  running  averages 

tLnp^tndSUKed11S0*°pe  ratlos<  We  have  also  correlated  the  1-sample  data  with 
temperature  by  least  squares  to  show  that  the  values  of  the  coefficients  are 

erforsra?e  ^Jnprn^nt  °n  ^  n'erhod  of  ayeraging,  although  the  associated 

trprnfe  !e«rCfS',.Sp?r:aP{rJer™oSSSr^f;0n  W,th  TEn9l«nd.  «• 

'  '•60]PPt/“C  for 
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TABLE  III  -  Comparison  of  measured  temperature  coefficients  with  those 


computed  assuming  formation  of  cellulose  at  equilibrium. 


COMPARISON  OF  MEASURED  and  COMPUTED*  TEMPERATURE 
COEFFICIENTS  OF  CARBON  AND  OXYGEN  IN  TREE  RINGS  IN 
PARTS  PER  THOUSAND  PER  DEGREE  CENTIGRADE. 


Computed* 

Assuming 

Equilibrium 

Measured 
in  this 
Paper** 

1 8 

for  exchange  of  0  between 
water  and  cellulose 

) 

0.36 

1 8 

for  exchange  of  0  between 
atmospheric  oxygen  and 
cellulose 

0.96  | 

3.7  +  0.5 

1 8 

for  exchange  of  0  between 

C02  and  cellulose 

1.14 

i 

for  exchange  of  C  between 

C02  and  cellulose 

0.36 

2.4  +  0.4 

for  exchange  of  D  between 
«ater  and  the  O-H  bond  in 
cellulose 

2.0  ' 

for  exchange  of  D  between 
water  and  the  C-H  bond  in 
cellulose 

0.4  . 

76  +  10 

*The  computed  values,  taken  from  Reference  1,  refer  to  isotope  exchange  at 
equilibrium  and  therefore  represent  lower  limits  of  the  real  values  of 
exchanges  which  may  depend  on  dynamic  effects.  The  measured  values  are 
larger  than  the  computed  values  and  thus  must  include  dynamic  effects, 
and  effects  from  formation  of  rain  water. 

**It  is  known  that  for  a  water  plant,  Euglena,  that  90 %  or  more  of  the 
oxygen  bound  in  cellulose  comes  from  water,  (Reference  12).  The  source 
of  oxygen  in  trees  has  not  been  determined,  but  probably  i  the  same. 
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i  10/H)  IN  PARTS  PER  THOUSAND 


Figure  4 


Plot  of  6d  vg.  6  in 


rain  and  snow  using  data  fron,  Ref 
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ABSTRACT 

We  consider  whether  the  long  tarn  perturbation  of  radiocarbon  dates, 
uhich  is  known  to  be  approximate! y  a  sin  function  of  period  about  B000 
years  and  amplitude  of  about  B*  peek-to-peek.  could  have  been  caused  in 
any  major  part  by  vulcanic.  Ue  conclude  that  this  ie  not  the  case. 

On  the  contrary,  present  day  volcanoes  are  a  far  less  important  source 
of  inert  C02  (about  100  fold  leas)  than  is  man's  burning  of  fossil  fuels 
ehich  haa  caused  the  Sussa  dilution  of  about  2%. 


I.  Introduction 

The  emission  of  C02  by  volcanoes  could,  in  principle,  reduce  the 
concentration  o T  radiocarbon  in  atmospheric  carbon  dioxide  and  thus 
cr educe  errors  in  the  radiocarbon  dates  for  the  past  millenie.  However 
me  corrections  are  opposite  in  sign  to  the  effects  that  are  substantiated 
-y  the  geologic  evidence.  Namely,  the  volcanic  explanation  mould  require 
that  over  the  laat  eight  millenie,  volcanic  activity  mould  have  had  to 
increase  substantially  up  to  the  present,  contrary  to  the  geologic 
evidence.  Considerable  evidence  about  the  nature  of  such  effects  is  et 
hand  from  the  etudiee  or  the  perturbations  of  the  radiocarbon  age  caused 
-/  the  burning  of  fossil  fuel  and  coneequent  amission  of  carbon  dioxide, 
wnich  like  volcanic  gasea,  contains  no  radiocarbon. 

II.  Material  Balance 

As  first  suggested  by  Suess1,  the  C02  generated  by  combustion  or 
fossil  fuels  in  the  period  1B5D-19S0  A.O.  mould  have  reduced  the  C14/^2 
ratio  in  the  atmoaphere  by  a  full  10*  had  ell  the  fossil  C02  remained  in 
the  atmosphere.  The  actual  observed  dilution  is  2.0  ♦  0.3*.  This  smaller 
effect  most  likely  ia  caused  by  the  reletively  short  residence  time  of 
C02  in  the  atmosphere  before  absorption  into  the  eee,  a  time  variously 
estimated  aa  being  batmeen  10  and  30  years1’2.  In  the  time  or  100  years, 
“Mich  is  much  longer  than  the  mixing  time,  most  of  the  foseil  CO  mill 
"eve  diaaolved  in  the  sea.  2 


D-2 


APPENDIX  E 

PRODUCTION  OF  RADIOCARBON  IN  TREE  RINGS  BY  LIGHTNING  BOLTS 


E-l 


VOL  78,  NO  26 


JOURNAL  OF  GEOPHYSICAL  RESEARCh 


SEPTEMBER  10,  1973 


Production  of  Radiocarbon  in  Tree  Rings  by  Lightning  Bolts 

L.  M.  Libby 

It  &  D  Associates,  Santa  Monica,  Calijornia  901,02 
School  oj  Engineering  ami  Applied  Science,  University  oj  Calijornia 
Los  Angeles,  Calijornia  90024 

H.  R.  Lukens 

C'ulJ  Radiation  Technology,  GulJ  General  Atomic 
San  Diego,  Calijornia  92112 

Lightning  holts  ujjpoar  to  hr  able  to  produce  nrutrons  by  processes  associated  with  the 
inci  Ural  ion  of  ions  to  an  amount  estimated  tc  8c  about  a  percent  of  the  cosmic  ray  neutron 
production.  Tin:  corresponding  effect  on  radiocai.  n  in  tree  rings  can  explain  the  well-known 
short-term  secular  errors  in  radiocarbon  dates  as  being  caused  by  climatological  fluctuations 
in  frequency  of  lightning  storms. 


The  short-term  secular  variations  of  radio¬ 
carbon  in  tree  rings  have  a  time  scale  of  about 
100  years  and  an  average  amplitude  of  a  few 
percent.  They  have  been  attributed  to  fluctua¬ 
tions  in  the  solar  cosmic  ray  intensity  [Sum, 
11iti5J  Tbe  likelihood  that  they  could  have  been 
caused  by  changes  in  the  volume  of  volcanic 
emissions  has  been  discounted  [Libby  and 
Libby,  1972| 

In  the  present  communication  we  discuss  “C 
production  in  thunderstorms  and  show  that  its 
present  day  magnitude  may  lie  a  few  percent  of 
"C  production  by  eosinie  rays,  so  that  hundred- 
sear  fluctuations  in  the  frequency  of  thunder¬ 
storms  rould  explain  the  short-term  secular 
v  a  mi  Hons  in  radiocarbon  in  tree  rings. 

Neutrons  have  Iwcn  observed  to  be  produced 
in  I  diorator)  electric  discharges;  therefore  it 
mav  he  l hat  they  are  produced  in  natural  dis¬ 
charges.  namelv,  lightning  bolts. 

In  a  lightning  storm,  as  described  by  Loeb 
I  1954 1  there  are  electric  fields  up  to  1500 
'oils  cm  giving  nse  to  lightning  bolts  10*  kw 
hours  in  magnitude,  with  average  electron  cnergv 
of  — •  150  kev  If  a  Boltzmann  distribution  is 
assumed,  a  finite  fraction  (10‘)  of  electrons 
have  energies  greater  than  2  2  Mov.  and  a  inueli 
smaller  fraction  (10  ,J)  of  electrons  have  ener¬ 
gies  greater  than  8  Mev  Furthermore,  other 
acceleration  mechanisms,  chiefly  pinch  effert, 
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seein  to  he  operating  and  thus  producing  copi¬ 
ous  quantities  of  neutrons  of  ~2.5  Mev,  ac¬ 
cording  to  experimental  observation. 

In  an  experiment  [Stephanakis  et  al.,  1972] 
in  which  electric  discharges  were  created  be¬ 
tween  electrodes  connected  with  nylon  threads, 
the  nylon  exploded;  electrons  and  heavy  ions, 
which  were  accelerated  to  approximately  equal 
energies,  on  the  average  about  10  kev,  and 
neutrons  of  about  2.5  Mev  were  produced. 
Nylon  threads  enriched  in  deuterium,  when 
they  were  similarly  exploded,  produced 
copious  2.5-Mcv  neutrons  from  the  d-d  reac¬ 
tions;  thus  it  was  shown  that  the  heavy  ions 
were  being  accelerated.  The  yield  of  2.5-Mev 
neutrons  from  exploded  hydrocarbon  threads  of 
natural  hydrogen  abundance  was  determined  as 
abi  ut  10*  neutrons  for  an  electric  discharge  of 
11  watt  hours,  in  a  discharge  a  few  centimeters 
long. 

For  a  lightning  bolt  of  10*  watt  hour  \Locb, 
1954,  p.  348]  the  number  of  neutrons  is  ex¬ 
pected  to  be  10".  Specifically,  this  production 
is  arrived  at  ns  follows. 

The  yield  V  of  neutrons  in  a  pinch  is  expected 
to  var\  as  the  square  of  the  density  N  of  tons 
and  as  the  volume  F  and  lifetime  t  of  the 
pinch,  so  that 

X±  =  (KlY  Xx 

F„  "  \NJ  l\.  r„ 

where  the  subscripts  l  and  v  designate  light- 
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ntng  bolts  and  exploding  wires,  respectively. 
We  take  ti  =  500  /. tsec ,  t„  =  50  nsec,  N,  «- 
10  *  g/cm’,  A.  ~  1  g/cm’,  V,  ~  10,T  cm*,  and 
V .  ~  10 1  ctti’,  from  which  it  follows  that 

y,/Y.  =  10*.  If  Y.  =  10*,  then  Y,  ~  10“ 

On  the  average  around  the  world,  R.  E. 
Holier  (private  communication,  1972)  has 
evaluated  the  total  number  of  lightning  bolts 
per  second  to  be  bet  ween  100  and  300,  which, 
if  Loeb’s  calibration  of  lightning  bolts  is  as¬ 
sumed,  would  lead  to  an  estimate  of  abollt 
1  X  10"  to  3  X  10"  neutrons  created  in  the 
atmosphere  per  second  over  an  area  of  5  X 
10"  cm*.  Holier  warns  that  this  calibration  is 
probably  too  high  by  at  least  a  factor  of  10. 
We  therefore  reduce  it  by  this  factor  and  arrive 
at  a  neutron  production  from  storms  of  (2-fl) 
X  10  *  neutrons/cin’  sec,  as  compared  with  the 
cosmic  ray  production  of  2.6  neutrons/cm*  sec, 
or  a  perturbation  of  about  1%. 

fin  ess  notes  that  a  comparison  of  recent  radio¬ 
carbon  fluctuations  in  tree  rings  with  eunsfiot 
numbers  since  1000  A.D.  suggests  that  a  corre¬ 
lation  exists  between  them.  It  is  not  unreason¬ 
able  that  frequency  of  thunderstorms  should  be 
correlated  to  sunspot  numbers  and  should 
modulate  the  production  of  atmospheric  radio¬ 
carbon  by  durtric  discharges;  thus  the  correla¬ 
tion  found  by  Suess  in  tree  rings  wuuld  be 
explained. 

It  is  barely  possible  that  an  enhancement  of 
radiocarbon  in  trees  on  a  mountar  ridge  fre- 
qinutlv  struck  by  lightning  could  be  measured. 
The  natural  rate  of  "C  decay  is  15  dpm  g  C, 
cor  re.- ponding  to  a  content  (taking  the  lifetime 
as  ,H265  years  or  4.1  x  10*  min)  of  6  x  10,u  "C 
atoni.a/g  (  .  Let  us  ask  for  a  1%  perturbation, 
namely,  that  t.  x  10*  “C  atoms/g  C  be  pro¬ 
duced  by  lightning;  in  dense  dry  (dead )  wood, 
this  corresponds  to  2  x  10*  "*C  atoms/g  wood. 

The  nitrogen  content  of  such  wood  is  about 
2  X  10  *  by  weight.  Neutrons  entering  wood 
are  captured  «•!  her  by  hydrogen  -j-  by  nitrogen. 
The  fraction  captured  hy  nitrogen  is  given  by 

t  =  ANffN/(  -4N<TN  +  AH(Tl|) 

where  ,4„  and  ,1„  are  the  atomic  fractions  of 
N  and  II  and  where  tr*  =  aN(«‘‘N  -*  p" C)  = 

1  7  x  10  "  cm*  and  cr„  =  a„(nH  -*  yD)  = 

0  3  x  10  *•  end 
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Taking  the  schematic  formula  for  wood  to  be 
CH,0  and  considering  capture  in  a  fallen  tree 
with  negligible  water  content,  such  as  Bristle- 
cone,  we  find  that  F  ~  0.02, 

To  create  1%  excess  of  "C  in  wood  then 
requites  capture  of  1010  neutrons/g  wood,  or 
10“  neutrons  in  a  tree  weighing  a  ton.  Thus, 
to  produce  this  amount,  several  hundred  light¬ 
ning  bolts  must  strike  nearby  trees  and  rocks. 
It  follows  that  fallen  trees  on  an  exposed  ridge 
where  electrical  disch.  rge  is  common  in  every 
storm  just  might  show  this  enhancement.  This 
"C  may  be  chemically  separable  from  cellulose 
and  lignin,  since  the  hot  atom  chemistry  of 
carbon  is  unlikely  to  introduce  the  "C  into 
these  molecules  without  shattering  them.  This 
possibility  could  be  checked  by  irradiating  wood 
with  neutrons. 

A  caution  should  be  added  that  the  foregoing 
estimates  of  "C  production  may  be  in  error  in 
'■ither  direction.  Namely,  the  value  of  the  aver¬ 
age  lightning  bolt  energy  [Loeb,  1954]  we  have 
used  although  it  was  decreased  by  a  factor  10 
may  still  be  too  high.  On  the  other  hand,  the 
numbers  of  neutrons  observed  in  experimental 
discharges  [SfepAanoAis  et  al.,  1972]  scatter 
over  a  range  of  up  to  5  X  10*  neutrons  of  2.5 
Mev,  and  there  may  be  much  larger  numbers  of 
lower-energy  neutrons  produced,  ell  of  which 
are  able  to  produce  "C  by  eapture  into  nitrogen. 
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Whon  recent  dnta  on  the  rate  of  addition  of  CO,  to  the  atmosphere  as  a  result  of  bacteriallv 
caused  decay  of  organic  carbon  produced  by  the  biosphere  are  used,  it  is  found  that  variations 
of  climate  acting  on  the  biosphere  can  cause  fluctuations  of  up  to  100  years  in  amplitude  in  the 
corrections  to  radiocarbon  dates. 


This  paper  discusses  the  effect  on  radiocarbon 
dates  of  variations  in  the  amount  of  carbon 
stored  globally  as  organic  material  that  origi¬ 
nated  in  the  biosphere.  Using  new  measure¬ 
ments  of  global  evolution  of  COt  from  the 
decay  of  stored  organic  carbon,  we  find  that 
tlio  short-term  ‘wiggles'  [Stress,  J970]  in  the 
amount  of  radiocarbon  in  tree  rings,  of  an  am¬ 
plitude  of  about  100  years  peak  to  peak,  can 
b«  explained.  The  present  discussion  updates 
and  quantifies  the  original  trc.’‘ment  of  this 
subject  by  Libby  [1952], 

Discussion 

lu  the  process  of  making  bio-organic  matter, 
plants  store  carbon  depleted  with  respect  to 
atmospheric  carbon  by  about  2.7 %  in  "C 
[Kroopnick  et  al„  1970;  Gordon  and  Williams, 
1970;  Maugh.  1972]  and  therefore  by  about 
in  UC  [Rafter,  1955],  If  the  amount  of 
depleted  carbon  stored  in  the  biosphere  in  the 
past  has  been  different  from  what  it  is  now 
radiocarbon  dates  for  the  past  will  reflect  an 
error.  The  size  of  error  is  estimated  in  this 
paper  and  is  found  to  be  about  100  years  maxi¬ 
mum. 

The  carbon  in  the  bio-organic  material  with¬ 
drawn  from  the  atmosphere  and  stored  m  the 
plant  for  the  life  of  the  plant  [ibis  an  additional 
time  of  storage  as  dead  organic  carbon  is  even¬ 
tually  released  by  bacterial  decay  and  returned 
to  the  atmosphere  mainly  as  CH„  In  the  atmo¬ 
sphere  methane  quickly  oxidizes  to  CO,  as  an 
end  product  At  present  the  global  quantity  of 
CO,  originating  from  plant  decay  and  returned 
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to  the  atmo.'tiherc  L  estimated  [Maugh,  1972] 
as  about  1.8  a  !<>"’  tons  CO,/yr,  or  1.0  X  10'# 
g  C/cm1  yr.  Let  the  total  time  of  withdrawal 
from  the  atmosphere  be  r  years  Then  it  follows 
that  the  amount  of  C  presently  stored  as  or¬ 
ganic  carbon  is  1.0  x  10'J  rg/cm’. 

The  storage  time  r  may  vary  depending  on 
climate  and  vegetation;  for  example,  the  mean 
age  has  been  found  to  range  from  10  to  1001) 
years  in  New  Zealand  soil  [O'Brien  and  Stout, 
1972J.  The  mass  of  organic  carbon  is  mainly  in 
the  oceans,  however.  The  average  concentration 
of  dissolved  organic  carbon  in  the  oceans  l- 
measured  f Menzcl  and  Hythcr,  1970;  Ilohn- 
Ilanscn,  1970]  as  about  500  /ig/1  or  0.5  pg  C/g 
sea  water.  When  the  average  depth  of  the  oceans 
is  taken  as  2S0O  meters,  it  follows  that  0.14  g 
organic  C/cm’  is  stored  as  dissolved  material 

Since  1  X  10''  g  C/cm*  is  returned  to  the 
atmosphere  per  year  [Maugh,  1972],  the  aver¬ 
age  storage  time  is  about  140  years. 

An  independent  estimate  of  the  storage  time 
comes  from  the  following  biological  argument. 
More  than  OS^  of  all  organic  carbon  in  the 
oceans  is  dissolved  organic  carbon,  at  about  500 
/'gd  [ Uohn-llanscn ,  1970],  The  total  quantity 
of  this  dissolved  organic  carbon  in  the  entire 
depth  of  the  ocean  represents  about  100  times 
the  average  annual  production  of  organic  carbon 
in  the  sea  surface.  The  annual  production  is 
known  quantitatively  to  just  be  able  to  satisfy 
the  metabolic  requirements  of  ti  e  grazing  popu¬ 
lation  [Mrnsel  ami  Ryther,  1970],  Thus  the 
total  organic  carbon  in  the  oceans  is  lhe  amount 
presently  produced  in  about  150  years  and  l 
in  steady  state,  being  eaten  and  re-eaten,  th> 
final  state  being  gases  produced  by  bacterial 
action.  The  chemical  estimate  of  140  years  and 
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(ho  bioloRir.il  estiinato  of  about  100  years  arc 
in  reasonable  aRreement, 

In  comparison,  Libby  [1952],  in  explaining 
the  fundamentals  of  radiocarbon  dating,  arrived 
at  an  inventory  of  0.33  r  organic  C/cm*.  His 
estimate,  made  25  years  ago,  was  remarkably 
close  to  the  above  estimate  of  0.14  g  C/cm * 
derived  from  new  data. 

The  total  atmospheric  CO,  is  about  2.33  X  10” 
tons,  or  0.13  g  C/cm’.  The  value  of  the  atmo¬ 
spheric  residence  time  before  absorption  into 
ihc  seas  has  recently  been  updated  [Gullikson 
and  Nydal,  1972]  to  about  8-12  years.  In  the 
absence  of  universal  agreement,  for  convenience, 
we  round  it  off  here  as  10  years.  Then  2.3  X  10" 
tons,  or  0  013  g  C/cm’,  are  exchanged  with  the 
oceans  each  year. 

In  past  times  the  total  mass  of  stored  organic 
carbon  may  have  been  larger  or  smaller  than 
it  is  now,  depending  on  the  climate.  Let  us 
define  as  the  norm  the  amount  of  carbon  pres¬ 
ently  stored  and  define  a  time  dependent  factor 
a  by  which  the  organic  carbon  reservoir  may  be 
increased  (<*  >  1),  as  in  a  lush,  tropically 
global  coal  age,  or  decreased  (a  <  1),  as  in  an 
ice  age.  We  assume  that  tbc  total  atmospheric 
CO,  has  always  remained  tnc  same. 

Since  the  organic  ea  bon  reservoir  is  about 
5 ^  depleted  in  radio, arbon,  the  atmospheric 
reservoir  is  correspondingly  enriched  by  about 
5%  multiplied  by  the  ratio  of  the  mass  of  stored 
organic  carbon  to  tiie  total  carbon  in  both  reser¬ 
voirs,  namely,  the  atmosphere  and  the  organic 
reservoir.  That  is,  the  enrichment  <  is  given  by 

*  =  0.05| 0. 14a/[0. 13  +  0.14(1  -  0.05)a)| 

‘  =  0.05(a/(0,95a  +  0.93)]  O 

where  «  is  the  ratio  of  the  mass  .1/  of  organic 
carbon  at  the  time  when  the  material  was  alive 
to  the  mass  M*  in  the  vear  1972  A  D  a  = 
M/M* 

The  enrichment  <  can  be  positive  or  negative 
with  respect  to  the  present.  Now  in  1972  A.D. 
we  define  a  =  1.0  and  <  =  ,*  =  0.02fifi.  The 
error  m  the  radiocarbon  date  witli  respect  to 
the  norm  is  given  by  t lie  enrichment  difference 
times  the  mean  radiocarbon  lifetime  of  8208 
years: 

err,"-  =  («  ~  **)  X  82G8  years  (2) 
The  age  correction  depends  on  the  numerical 


value  of  5%  depletion  of  "C  in  plants  in  com¬ 
parison  with  JO,  in  the  atmosphere.  This  quan¬ 
tity  in  turn  depends  on  the  temperature  T  at 
which  the  bio-organic  material  grew  and  itself 
has  a  temperature  coefficient.  The  coefficient 
for  ”C/"C  in  marine  plankton  has  been  mea¬ 
sured  as  0.35  ppt/°C  [Be  et  al ,  1965]  and 
independently  as  0.5  ppt/°C  [Eadie,  1972],  It 
may  be  the  same  or  larger  for  the  average  of 
■*ll  ptants  [Libby,  1972],  In  the  absence  of 
more  measurements  we  assume  it  to  be  the 
larger  of  the  two  experimental  values,  namely 
0.50  ppt/°C,  because  it  seems  nriikely  that  the 
true  value  is  exactly  equal  tr  the  minimum 
value  allowed  thermodynamically  [Libby,  1972], 
which  is  0.35  ppt/°C. 

Then,  on  the  basis  of  the  foregoing  assump¬ 
tions,  the  temperature  coefficient  will  be  1.0 
PPt/  C-i  so  that,  if  the  climate  were  much 
warmer,  say,  by  10°C,  the  "C  depletion  of  the 
stored  organic  carbon  would  decrease  from 
0.050  to  0.040,  corresponding  to  an  age  correc¬ 
tion  of  (0.050  -  0.040)  x  8208  years,  or  83 
years  too  old.  Similarly,  if  (he  temperature  fell 
10°,  the  depletion  would  increase  to  0.060, 
corresponding  to  an  age  correction  of  83  years 
too  young,  even  without  any  change  in  the  mass 
of  the  biosphere,  the  dependence  on  which  is 
set  forth  in  (2). 

Rut  still  another  effect  enters:  that  caused  by 
the  presenee  of  ocean  carbonate.  The  sea  con¬ 
tains  dissolved  carbonate  that  is  fractionated 
in  "C/”C  in  comparison  with  that  in  the  atmo¬ 
sphere,  the  amount,  of  dissolved  carbonate 
depending  on  the  temperature  at  the  sea-air  inter¬ 
face.  Sea  carbonate  is  enriched  in  "C  with  re¬ 
spect  to  atmospheric  CO,  and  is  less  so  as  the 
temperature  increases,  so  that  the  result  is  to 
provide  a  small  correction  that  partly  cancels 
the  effect  of  fractionation  in  stored  organic 
carbon. 

This  port urbai ion  conics  only  from  carbonate 
in  the  ~  IdO-met  er-thick  surface  mixed  layer 
of  the  ocean  thai  contains  [Craig  et  al.,  1972) 
about  0.15  g  C'/cni3,  which  is  approximately 
equal  to  that  in  the  atmosphere  and  to  that 
m  the  organic  carbon  dissolved  m  the  total 
depth  of  the  sea.  In  1  lie  ease  of  organic  carbon 
the  total  depth  iff  the  sea  is  involved  because 
bacterial  decomposition  occurs  at  all  depths, 
producing  (  II.  and  <’(),  with  which  the  sea 
water  is  saturated,  so  that  these  bacterial 
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gases  are  bubbling  up  from  all  depths.  In  com¬ 
parison,  the  sea  is  not  saturated  with  CO.  at 
any  depth  [Craig  et  ai,  1972], 

Carbon  isotope  fractionation  by  CO,  absorp¬ 
tion  at  the  air-water  interface  has  been  mea¬ 
sured  [Degens  and  Deuser,  1967],  and  the  ratio 
“C/“C  lias  been  found  to  vary  from  9.2  to  6.S 
ppt  over  the  temperature  range  0°-30°C.  The 
corresponding  fractionations  of  “C/’C  versus 
temperature  are  shown  in  Table  1  and  are 
obtained  by  multiplying  the  fractionation  of 
”C/”C  by  2. 

The  corrected  enrichments  (or  depletions)  of 
"C/”C  in  the  atmosphero  for  climates  10°C 
warmer  (or  colder)  than  those  at  present  arc 
given  in  Table  2,  where  it  is  assumed  that  the 
average  global  temperature  at  present  is  20°C. 
The  corrections  caused  by  ocean  carbonate  are 
shown  to  lie  about  20%  of  the  temperature 
variations  of  l‘C/,*C  in  organic  carbon  and  of 
opposite  sign. 

The.  equation  including  both  corrections  ean 
now  be  written 

t  =  na/[(l  -  n)a  +  0.93]  ^ 

n  =  [0.05  -  8.5  X  10'-(T  -  T*)] 

The  computed  corrections  are  shown  in  Table  2 
to  he  of  the  order  of  ±50  years,  namely,  of  a 
total  amplitude  cf  about  100  years. 

We  conclude  that  the  experimentally  deter¬ 
mined  fluctuations  in  the  correction  curve  for 
radiocarbon  dales  [.S' ness,  1970]  as  deicmuncd 
in  tree  rings  may  be  caused  by  changes  in  the 
climate  that  perturb  the  storage  mechanism 
discussed  here.  A  short  period  of  climate  either 

TABLE  1.  Corrections  Caused  by  Fractionation  of 
C  in  Surface  Ocean  Carbonate  and  in  Dissolved 
Organic  Carbon  Versus  Average  Temperature 
of  the  Sea  Surface,  a 


Sea  Carbonate 
in 

Mixrj  Layer 

Dissolved  Organic 
Carbon  in  Sea 

Temperature  t 

*C 

1  ♦  il.018(i* 

1  -  0.070 

0 

1  ♦  0.017 

1  -  0.060 

10 

1  •  U.(J15t 

1  -  0.050 

20 

1  *  0.0138 

1  -  0.040 

30 

“  is  the  ratio  ot  (|I*C/1JC)  in  aqueous  car- 
boiiate  to  (,I-C7,JC)  in  atmospheric  CO,. 

•As  given  by  Craig  et  al.  11972) .  2 
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Corrections  to  Radiocarbon  Dates  in 
Bio-Organic  Materials  as  a  Function  of 
Average  Temperature  and  Biomass 


a  ■  H/H' 

Correction,  years 

10*C 

20*C 

30*C 

0.5 

-48 

-7J 

-111 

1.0 

38 

0 

-40 

2.0 

127 

73 

18 

296 

214 

144 

Corrections  assume  that  the  changed  climate 
endures  for  100  years  or  more  and  include  the 
corrections  caused  by  ocean  carbonate  in  the 
mixed  layer  of  the  sea. 

warmer  or  colder  than  that  at  present  could 
cause  the  atmospheric  radiocarbon  to  differ 
from  the  norm,  with  a  corresponding  error,  posi¬ 
tive  or  negative,  in  radiocarbon  dates  of  bio- 
organic  material  for  those  periods.  By  short  is 
meant,  periods  of  time  about  as  long  as  the 
average  storage  time  of  about  100  years.  This 
storage  time  itself  might  have  changed  with 
tune  if  the  oceans  had  changed,  but  the  tem¬ 
perature  of  the  sea  below  the  thennocline  is 
only  slightly  above  4°C  and  the  time  of  over¬ 
turn  is  about  1000  years,  so  that  no  significant 
change  could  have  occurred  in  the  last  8000 
years,  the  period  for  which  the  fluctuations  in 
tree  ring  radiocarbon  have  been  mea-nred. 

Support  for  the  hypothesis  of  a  fluetuaiiou  m 
the  amount  of  stored  organic  carbon  in  the  past 
is  suggested  by  measurements  on  CO  m  the 
flreenland  and  Antarctic  icecaps  [Cavanagh  cl 
1972J.  Specifically,  these  measurements  find 
concentrations  of  CO  as  high  as  0.3  ppm  in  air 
trapped  in  the  ice  lor  ice  samples  up  to  2500 
years  old,  whereas  modern  CO  concent  rat  ions 
in  clean  air,  remote  from  cities,  a*v  0.04-0  20 
ppm  at  ground  levels  |  Jajje,  1972 1  and  0.0S 
ppm  at  a  4-kni  altitude  and  derrea.-e  to  0.04 
ppm  at  it  15-kill  altitude  f  Bradford,  19721  Since 
the  largest  global  source  of  CO  is  oxidized 
methane,  these  measurements  suggest  that  the 
reservoir  of  stored  bio-organic  material  has 
varied  in  the  past. 
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Tho  principle  of  isotope  thermometry  is  that  isotope  ratios  are  temperature  dependent. 
It  is  proposed  that  isotopic  ratios  in  bio-organic  material  may  have  independent  temperature 
coefficients,  so  that,  by  measuring  ratios  for  several  elements,  e.g.,  hydrogen,  carbon,  and 
oxygen,  it  may  be  possible  to  show  that  a  temperature  change  occurred.  Assuming  that 
cellulose  is  formed  almost  at  equilibrium,  temperature  coefficients  for  C,  H,  and  O  are  com¬ 
puted  and  found  to  be  of  measurable  size.  The  coefficients  computed  for  ”C/‘’C  and  D/H 
agree  with  measured  values.  Thus  indications  are  that  temperature  changes  in  past  climates 
can  be  measured  in  old  tree  rings  by  using  this  principle  of  multiple  thermometers. 


The  principle  of  isotope  thermometry  is  that 
the  isotopic  ratios  are  temperature  dependent. 
There  are  two  kinds  of  temperature  depen- 
'ence,  kinetic  and  equilibrium.  The  theory  of 
netie  isotope  separation  is  well  established, 
or  example,  in  the  gaseous  diffusion  process  for 
separation  of  the  uranium  isotopes,  separation 
is  caused  by  a  difference  in  the  rates  with  which 
the  two  species  pass  through  holes  in  a  barrier. 

In  case  of  equilibrium,  separation  is  caused  by  a 
difference  of  difference.',  namely  the  difference 
between  separation  by  the  forward  rate  and 
separation  by  the  reverse  rate,  the  two  rates 
being  equal.  Distillation  of  liquid  in  a  closed 
container  is  an  example  of  separation  at  equi¬ 
librium,  while  vaporization  into  a  vacuum  is 
kinetic. 

In  distillation  of  a  liquid  containing  two  or 
more  elements  in  t lie  molecule,  the  correspond¬ 
ing  isotope  thermometers  are  functionally  re¬ 
lated.  because  the  corresponding  vapor-pressure 
ratios,  which  determine  the  separations  on  evap¬ 
oration  and  condensation,  depend  on  the  same 
force  constants  tying  them  together. 

In  bio-organic  material,  however,  there  may 
he  ai  least  three  Useful  isotope  thermometers, 
namely  I)/ll,  “C/'T,  and  '■0/'"0.  The  basic 
reunion  is  pholo-ynthesis  wiih  production  of 
cellulose,  here  written  schematically  as  a  basic 
module, 
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H  -  C  -  OH  or  CHaO 

I 

according  to 

COa  +  HjO  ->  CHaO  +  Oa 

In  this  reaction,  all  three  elements  may  yield 
isotope  thermometers,  and  they  may  be  inde¬ 
pendent  of  each  other;  that  is,  the  ratio  ‘"O/'O 
is  no  longer  functionally  related  to  D/H,  be¬ 
cause  diatomic  oxygen  is  being  evolved,  and  so 
rotational  and  vibrational  energies  of  the  0-0 
bond  are  involved,  and  because,  in  cellulose, 
oxygen  is  bound  10  carbon  as  well  as  hydrogen. 
For  the  same  reason,  the  temperature  coefficient 
of  “0/"0  and  nC/"C  may  not  be  functionally 
related.  Finally,  the  ratio  ”C/”C  should  be  in¬ 
dependent  of  D/H,  because  carbon  and  hydro¬ 
gen  are  initially  in  different  molecules  with  dif¬ 
ferent  vibrational  frequencies. 

If  variations  are  measured  for  only  one  iso¬ 
tope  ratio,  the  rclalion  to  paleotemperature  is 
not  proved.  However,  if  ratios  are  measured  for 
two,  three,  or  more  independent  isotopes  in  the 
same  organic  material,  it  is  possible  that  the  re¬ 
lation  to  temperature  can  he  firmly  established 
and  quantified,  and  this  is  the  principle  point  of 
the  present  discussion. 

Over  the  small  temperature  interval  (0°C  < 
t  <  +  30°C)  within  which  ehmate  could  have 
varied  without  killing  terrestrial  life  as  we  know 
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it,  t he  dependences  of  isotope  ratios  on  tem¬ 
perature  ean  be  assumed  to  be  linear  for  bio- 
orgame  material.  That  is,  for  isotojie  ratios  S', 
where  j  =  2,  13,  15,  18,  34,  41,  ••• ,  referring 
to  D,  "C,  UN,  **0,  “S,  "K,  •••  ,  the  tempera¬ 
ture  dependence  ean  be  written  as 

S'  =  a't  +  b> 

(a i  and  b'  constants) 

where  a’  is  the  temperature  coefficient  and  t  is 
temperature.  Then  a  change  A'  is  related  to  a 
change  AS'  by  (for  a  given  compound) 

A(  AS1  AS13  AS14  AS18 
a  a  a  a 

The  mult ijile  overdetermination  of  temperature 
by  using  a  set  of  many  thermometers  may  make 
it  possible  to  show  with  some  level  of  confidence 
that  a  temperature  change  did  occur.  It  is  also 
possible  that  new  effects,  not  caused  by  tem¬ 
perature,  will  he  demonstrated  by  the  method  of 
multiple  thermometry,  if  some  of  the  above 
ratios,  but  not  all  of  them,  show  equality'.  Prob¬ 
ably  the  specific  relations  to  temperature  will 
have  to  be  calibrated  at  each  geographical  posi¬ 
tion  to  take  into  account  local  variations  in 
water  and  CO,. 

Finally,  bio-organic  data  bases  can  be  radio¬ 
carbon  dated,  so  that,  if  evidence  for  tempera¬ 
ture  changes  is  stored  in  them,  the  dates  when 
the  changes  occurred  ean  he  evaluated. 

Experimental  Evidence  for  Temperature 
Dependence  of  Isotope  Patios  in 
Dio-Oruanic  Material 

Although  animal  hio-organie  material  is  not 
directly  derived  from  photosynthesis,  neverthe¬ 
less  animals  are  largely  what  they  eat,  so  that 
bio-organic  material  from  both  plants  and  ani¬ 
mals  may  be  useful  in  isotope  thermometry. 
Some  of  the  evidence  [Craig,  1953]  on  isotope 
thermometers  in  bio-organic  material  is  as 
follows. 

The  "C/’V  ratio  iu  eombustible  organic  mat¬ 
ter  in  cores  from  sea  bottoms  shows  variations 
of  as  much  as  0  V,0  with  depth.  These  have 
been  measured  and  attributed  to  temperature 
variation  in  the  sea  surface  with  time  by  Rogers 
iiml  Kuoit.s  (I960].  Also,  in  a  study  of  marine 
plankton  [ Sacked  et  at.,  1905]  the  UC  concen¬ 
tration  was  found  to  be  depleted  by  0  per  mil 


where  surface  waters  are  near  0°C  relative  to 
samples  collected  where  surface  temperatures 
are  about  25°C.  These  observations  on  both 
plant  and  animal  plankton  suggest  that  the 
temperature  coefficient  of  UC/“C  in  bio-organic 
material  is  ~0.24%r/oC  in  the  temperature 
range  0°-25°C. 

In  further  work  [Degens  el  at.,  1968],  the 
temperature  coefficient  of  "C/“C  in  marine 
phytoplankton  grown  in  laboratory  cultures  was 
measured  as  ~0.35"r/°C  between  10°  and  30°C. 

Although  shell  carbonate  is  not,  strictly  speak¬ 
ing,  bio-organic  material,  nevertheless,  for  com¬ 
pleteness  we  note  here  that  variations  in  the 
‘*0/'*0  ratio  in  Fora  min  ifera  shells  from  Carib¬ 
bean  Sea  cores  have  been  measured  by  several 
workers  and  interpreted  in  terms  of  paleotem- 
perature  variations  [ Emihnni,  1966].  This,  in 
fact,  was  the  first  isotope  thermometer  devel¬ 
oped  by  Urey  [1947], 

It  may  be  that  other  i-otope  ratios  in  bio- 
organie  material  are  temperature  dependent, 
such  as  'W/'N  and  “S/^'S,  perhaps  increasing 
the  number  of  independent  thermometers  in  bio- 
organic  data  bases  to  five  or  six.  This  possibility 
should  be  investigated. 

The  Question  of  Equiluirium 

An  attempt  is  made  here  to  estimate  the 
magnitudes  of  the  temperature  coefficients  of 
the  several  isotope  thermometers  in  bio-organic 
data  bases.  This  can  be  done  using  thermody¬ 
namic  considerations,  if  one  assumes  that  the 
organic  material  is  manufactured  in  equilibrium 
with  its  surroundings.  This  is  of  course  a  major 
assumption,  but  one  that  does  not  seem  too 
unlikely.  For  example,  in  experiments  in  the 
laboratory,  the  isotopic  composition  of  CaCOs 
slowly  formed  from  aqueous  solution  was  found 
to  be  the  same  as  that  in  the  shells  produced  by 
aquatic  organisms  at  the  same  temperature 
[McL’rea,  1950] .  In  any  case,  this  assumption 
can  be  tested  for  plants  and  animals  grown 
under  laboratory  conditions. 

In  the  wild  state,  in  coastal  waters  and 
estuaries,  the  carbon  and  oxygen  isotopic  com¬ 
position  of  shell  carbonate  has  been  shown  to  be 
in  isotopic  equilibrium  with  bicarbonate  dis¬ 
solved  in  the  water  [Muni:  anil  Vogel.  lflCs], 
The  life  span  of  die  shell  fi-h  is  aboui  one  year, 
so  it  would  seem  that  equilibration  in  K‘lf- 
manufactiire  might  take  place  rapidly,  say  in 
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less  than  or  about  a  day.  Equilibration  between 
CO,  pas  and  sea  water  is  known  to  occur  in 
less  than  48  hours  (Li  and  Tsui,  1971],  and 
among  CO,  gas,  water,  and  precipitated  CaCO, 
in  less  than  G  days  [Dottinga  and  Craig,  19G9; 
Ernrirh  et  al.,  1970],  In  living  animals,  tritiated 
water  is  known  to  equilibrate  with  blood  serum 
of  rats  in  0  days  or  less  [ Gleasun  and  Freedman, 
1970],  Tims  the  evidence,  although  meager,  is 
that  equilibration  takes  place  rapidly  compared 
with  the  relevant  life  spans. 

In  shells,  the  carbon  and  oxygen  are  not  inde¬ 
pendent  thermometers,  because  they  are  tied  by 
the  same  vibrations  in  the  C-0  bond.  However, 
there  is  organic  material  in  the  shell  matrix  that 
could  provide  a  set  of  multiple  thermometers, 
and  also  could  be  radiocarbon  dated. 

The  question  of  the  degree  to  which  equilib¬ 
rium  obtains  in  formation  of  cellulose  is  interest¬ 
ing  in  itself.  If  equilibrium  seems  to  apply,  t be 
assumption  could  be  extended  to  computations 
of  production  of  other  substances.  If  equilibrium 
does  not  obtain,  and  instead  the  processes  of 
isotope  fractionation  are  kinetic,  then  the  tem¬ 
perature  coeflicients  should  he  even  larger  than 
those  estimated  by  assuming  equilibrium. 

Calculation  of  Partition  Function  for 
Cellulose 

For  a  reaction  in  which  there  is  isotope  ex¬ 
change.  the  equilibrium  constant  A',,  can  be  ex¬ 
pressed  m  terms  of  the  total  partition  functions 
Q*  and  Q  for  molecules  containing  the  heavy 
and  light  isotopes,  respectively,  as  follows.  We 
define  Q  as  the  product  of  rotational,  vibra¬ 
tional.  and  traii'lational  parts,  neglecting  vibrn- 
t ional-rot at ional  coupling : 

Q  =  Q,..<UQ.,  (i) 

For  solids  and  liquids,  Q„  =  1  and  Q,lb  is 
given  by 

Q*.n(T)  =  Z  I  J,  exp  v,n  ~j 

exp  (-St, /AT)  (n  =  0,  I,  2,  ••  •)  (2) 
which  can  be  written  as 

Qy,i(T)  =  II, g, [I  -  exp  (  —  hv,/kT)\~' 

•exp  (-At, /k T)  (3) 

We  shall  be  computing  partition  functions  at 


TABLE  1.  Stretching  Vibrations 
Assumed  for  Cellulose 


Bond 

Stretching 
Vibration , 
cm- 1 

1 

—  C  —  II 

t 

t  1 

2960 

-c  -  c  - 

1  1 

900 

-0  -  II 

3680 

-C  -  (OH) 

(1200)* 

O 

lO  1 

T 

O  i 

sH 

1580 

B  =  1.44560 

^Calculated  from  a  force 
constant  taken  as  1/2  of  force 
constant  for  C  =  0  [Herzberg , 
1945,  p.  195,  Tables  51  and  89]. 

— 300® K,  so  that  the  harmonic  oscillator  ap¬ 
proximation  is  reasonably  good.  Here  i=  the 
Boltzmann  con-iam,  t  is  the  absolute  tempera¬ 
ture.  v,  is  the  frequency  of  the  ;ih  appropriate 
vibration  of  degeneracy  g,.  n  is  i he  vibrational 
quantum  number,  and  At,  is  the  difference  in 
zero-point  energy  between  the  isotopes.  In  cellu¬ 
lose,  the  barriers  preventing  free  internal  rota¬ 
tion  tire  high,  so  that  rotations  deteriorate  into 
torsional  oscillations  or  bending  vibrations.  In 
the  present  treatment,  these  will  be  included  in 
Q hi.  as  part  of  the  vibrational  partition  function. 

The  stretching  and  bending  vibrations  appro¬ 
priate  to  eellulo.-e,  which  i<  a  chain  molecule  of 
module  (II-C-OI1).  are  listed  in  Tables-  1  and 
2.  Their  frequencies  differ  slightly  for  molecules 
containing  difTereni  isotopes,  because  the  iso¬ 
topic  mas-  enter.-  into  the  reduced  mass  /i  of 
tilt'  appropriate  oscillator  from  which  the  par¬ 
ticular  frequency  i-  computed,  according  to: 


where  K  is  the  force  constant  of  the  particular 
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viliration.  Substituting  in  (’.l)  the  frequencies 
from  Table  1,  anil  noting  that  the  torsional  fre¬ 
quency,  tfiO  cm  is  threefold  degenerate,  one 
computes  y(272°K)  =  1.2d()2  and  Q(2!).S°K )  = 
1.29M  for  the  module  {11-C'-()II}.  In  principle, 
the  I )- 1 1  isotope  elTeet  in  cellulose  involves  two 
possible  substitutions:  namely,  substitution  of  a 

TABLE  2.  Bending  Vibrations 
Assumed  for  Cellulose 


Bending 

Vibrat ion 

Bending 
Frequency , 
cm- 1 

A 

R—  C  —  Oil 

1 

920 

1 

R 

11  A 

1  V 

R —  C  — Oil 

I 

920 

1 

R 

II 

| 

R — C  —Oil 

U 

920 

ii 

1  K! - Si 

R —  C  —0  —  11 

1 

(700)* 

1 

R 

II 

1 

1 

R —  0  —  Oil 

\  l 

K 

(375 J  , 

‘/;V r.;:rrj  [1945,  p.  196,  Tables 
5 1  and  89  j . 

■'rilerzberg  ( 1945,  Tables  118  and 
1 19  J  . 


TABI.I.  3.  Pcrtitiun  Function:;  CcnputsJ 

I 

for  the  Module  of  Cellulose  (11 — (*  —  01 1 ) 

I 

and  Its  Isotopic  Modifications 


Part  it  ion 
lunct i on 

273°K 

2'JS"8 

<3d  iron) 

1.231)2 

1  .2‘J8 1 

c*(ncoh) 

■10,850 

52. oas 

C*  (IKJ01»  > 

3T>  .011 

2R.0M> 

<3 * ( '  3CII-0) 

i 

i  :.;n 

<3* (<ll.  I0(l) 

i  .  jfijo 

i .( '.:>•> 

deuieroii  for  the  hydrogen  linked  to  the  earbon 
atom  and  substitution  of  the  hydrogen  attached 
to  oxygen  In  llte  ca-e  of  11  and  1C  a'taelied  to 
thi’  pivotal  carbon,  the  appropriaie  reduced 
ina->  is  insensitive  to  substitution  of  II  by  D 
because  the  masses  of  11  and  1C  are  almost  in¬ 
finitely  large  relative  m  II  or  0.  The  partition 
functions  computed  for  the  various  i-otopie  sub¬ 
stitutions  of  the  cellulose  module  are  given  in 
Table  :t. 

Heeaiisi'  of  the  cni'icncss  inherent  in  assuming 
bond  frequencies  for  cellulose  to  be  equal  to 
tlio-c  known  in  rather  small  organic  molecule-, 
we  neglect  inharmonic  corrections'  and  compli¬ 
cations  such  a-  hydration  in  solution  or  effect 
of  hydrogen  bond  formation. 

The  fractionation  ratios  can  be  written  as 
follows.  For  the  oxygen  isotopes,  there  are  three 
equilibriums  corresponding  to  isotopic  exchange 
between  cellulose  and  each  of  the  ihrcc  oxygen- 
containing  molecules  involved  in  the  reaction: 


ID 


K(02)  = 


"K(C0,)  = 


,BK(IIaO)  = 


V(Cl\,"  0)  qC  (),) 

4)(CH2'*0)  0*(,lO,iO) 

Q*(C1I,'\))  Q(C'\h) 
q{cu2'°o)  (A'c'Vo) 

(A£n?’“°J  CHijAb 

<?((TIj,r'())  A!l2‘"0) 


(5a) 


(56) 

(5e) 


where  ()*  is  the  |iartitiou  function  for  the 
licavier  isotope.  For  the  hydrogen  isotopes,  two 
rqmhhrunus  enter  so  that 


KflK'OD) 

K(I)COH) 


0>*(iH'<n>)  Q(n,<» 

V«'HM  Q'(Ul)O)  fri 
(b> 

0*(I)COII)  Q(1I20) 

QtCUjU)  Q*(l  IDO) 
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and,  for  the  carbon  isotopes,  only  one  equilib¬ 
rium  enters,  that  with  CO,,  so  that 


"K  =  Q*(13ch2o)  Q(‘acoa)  ... 

<?(,aCHaO)  Q*(“COj)  (  J 

but  what  the  chemist  actually  measures  [Urey, 
1947]  is  an  isotopic  fractionation  a.  For  exam¬ 
ple,  for  oxygen  the  fractionation  factor  corre¬ 
sponding  to  equation  5a  is: 

{ 2[laOa1  +  [“0“Q11  [CHa180] 

{[“O'°0J  •+•  2[,0Oa]}  [CHa^O]  °  {8) 

Considering  the  reaction 

'“Oj  +  ‘"Oj  J=?  2["0“0]  (9) 

the  equilibrium  constant  is  4,  neglecting  a  tiny 
correction  for  isotope  effect  [Urey,  1947]: 


4  =  r0“0]V["0a][,‘0a]  (10) 

because  ”0,  and  “0,,  being  symmetrical,  have 
only  hal  as  ma;  y  rotational  states  as  the  asym¬ 
metrical  molecule  w0'"0.  Substituting  (10)  into 
(8),  the  fractionation  factor  reduces  to 


a‘8(0a) 


<?roar  <?*[CHa‘*o] 

<2*[‘"Oa],/2  <3lCHa“0] 


(ID 


Corresponding,  CO,  and  C"0„  being  sym'- 
metric,  have  only  half  the  rotational  states  as 
C“0'*0,  so  that 


o“(C0j) 


QIC'W3  Q*l CHal8Q] 
<?*[C**03],/2  0(CHa1(,0] 


(12) 


but  for  a(HCOD)  and  o(DCOH)  the  isotope 


TABLE  4  Partition  Function  Ratios 
Computed  for  Deuterated  Cellulose 


Partition  Function 

273*K 

298'K 

U*|l'COII]/C|IICOHj<4 

CMHCODj/ejllCOIlj 

33.206 

25.187 

28.462 

21.636 

K*|l'?Olsai/0|U,O|  j'J- 

16.503 

12.513 

32.7400 

24.9460 

The  partition  function  ratio  for  liquid  is 
obtained  by  multiplying  the  partition  function 
ratio  for  gas  by  the  0/11  ratio  in  liquid  divided 
by  that  for  gas.  It  is  talen  here  as  1.104  at 
273*K,  and  as  1.074  at  298*K  IMerZivat  et  al. 
1903) . 
oTable  3. 
bUre y  (1947). 


TABLE  5.  Partition  Function  Ratios  Computed 
for  Cellulose  Containing  Heavy  Carbon  13 


Partition  Function 

273*K 

298*K 

C* ( 1 3CH20) /J( 1 JCH20)  “ 

1.1575 

1.1435 

In  C*[13C02)/U(>2C02]  (gas)* 

0.19732 

0.17558 

C* I 1 3C02 1 /c( 1 2C02 )  (gas) 

1.2181 

1.1919 

“Table  3. 
bBottinga  (1967). 


effect  is  large,  so  that 


o(HCOD) 

=  Q»[HC0D1  /[2Q[Ha0]  +  QfHDOl) 
Q[CHaO]  \[Q[HD0]  +  2Q[DaO]j  . 

o(DCOH) 

=  Q*[DCOH]  /2Q[HaO]  +  QfHDOl) 
Q[CHaO]  \0iHDO]  +  2Q[DaO]J 


Here,  while  H,0  and  D,0  are  symmetric,  none 
of  the  deuterated  or  hydrogenated  cellulose 
modules  are  Finally,  because  "CO,  and  “CO, 
have  the  same  number  of  rotational  states,  also 
H,*0  and  H,“0,  the  fractionation  factors  are 


o'*(HaO) 


Q(Ha160)  Q*(CHa‘"0) 
Q*(Ha180)  Q(CHa,80) 

Q(l8COa)  Q*(13CHaO) 
Q*Cc02)  Q(,2CHaO) 


(14 

(15) 


For  the  isotope  ratios  D/H,  »C/"C,  and 
lt0/l*0,  we  can  now  compute  the  partition 
function  ratios  (Q*/Q)cu, o.  The  partition  func¬ 
tion  ratios  for  (Q*/Q) OI„,a  are  taken  from  Urey 
[1947],  those  for  (H,0)iiq  are  taken  from  com¬ 
putations  of  Bottinga  [1967]  and  those  for  (CO,),., 
from  Bottinga  [1967]  and  for  II  DO  and  d',0 
from  Urey  [1947]. 

The  ratios  appropriate  for  each  isotope  ex¬ 
change  are  listed  in  Tables  4,  5,  and  6.  By  sub¬ 
stituting  these  in  equations  (11)  through  (15), 
tho  fractionation  factors  have  been  calculated 
for  27d°K  and  29S°K.  The  corresponding  tem¬ 
perature  coefficients  are  h-ted  in  Table  7. 

The  temperature  coefficient  so  calculated  for 
"C/“C  of  0.36%o/°C,  agrees  with  the  coefficient 
measured  for  ”C/“C  in  plant  and  animal  ma¬ 
terial  by  Sackett  and  his  co-workers  [Sackett, 
et  al.,  1965;  Degens  et  al.,  1968],  so  that  the 
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TABLE  6.  Partition  Function  Ratios  Computed 
for  Cellulose  Containing  Heavy  Oxygen 


Partition  Function 

273°K 

298°K 

Source 

0*[CH218O]/<2[C11216O] 

1.2697 

1.2833 

This  paper 

t«*(18o2)/«?c16o2)  ]gas,s 

1.0923 

1.0818 

Urey  (1947] 

b  In  |«*[C1802]gas/(J[C1602]gas| 
|«*[c18o2]/e[c16o2]}  gasb 

0.  12530 

1.1336 

0.11108 

1.1175 

Bottinga  [1967] 

In  ^*[ll2180]gas/«|L,160]gas| 

(«*[h218o]/<j[h218o]}  gas 

0.06822 

1.0706 

0.06164 

1.0635 

Bottinga  [1967]a 

aThe  partition  function  ratio  for  liquid  water  is  obtained  by 
multiplying  the  partition  function  ratio  for  water  vapor  with  the 
180/160  ratio  in  liquid  water  divided  by  the  ratio  in  the  vapor.  This 
correction  is  taken  as  1.01150  at  273®  and  as  1.00930  at  298°K  [ Bottinga . 
1967,  p.  806]. 


assumption  of  equilibrium  appears  to  be  more 
or  less  valid,  and  one  can  hope  that  the  like 
calculations  for  oxygen  and  hydrogen  may  be 
meaningful. 

At  equilibrium,  the  temperature  coefficient  of 
HDO/H,0  is  312£e/®C  [Urey,  1947],  much 
larger  than  that  of  the  (D/H)  ratio  calculated 
for  cellulose,  so  that  measuring  it  in  a  geographi¬ 
cal  distribution  of  the  same  kind  of  trees  would 
measure  the  distribution  in  D/H  ratios  of  rain 
water,  and  measuring  it  in  a  chronological  se¬ 
quence  would  measure  time  dependence  of  the 
rain  water. 

Our  naive  assumptions  are  mainly  as  follows. 
We  have  assumed  that  cellulose  is  formed  in 
solution,  that  in  cellulose  rotations  have  degen¬ 
erated  into  torsional  vibration,  that  its  vibra¬ 
tions  are  like  those  in  small  organic  molecules, 
and  that  manufacture  of  cellulose  in  plants  is  an 
almost  equilibrium  process.  (  A  further  assump¬ 
tion  has  been  made  in  the  neglect  of  acoustical 
vibrations  of  cellulose  in  the  lattice  of  solid 
wood.  While  these  might  make  a  contribution 
to  the  fractionation  factors,  their  temperature 
coefficients  appear  to  be  negligibly  small,  judg¬ 
ing  from  the  estimates  of  Bottviga  [1907]  for 
solid  calcite;  he  points  out  that  the  lattice  con¬ 
tribution  to  the  temperature  eocllicient  is  small. 
Our  main  concern  in  the  present  [taper  is  with 
the  magnitudes  of  the  temperature  coefficients, 


rather  than  with  the  fractionation  factors  them¬ 
selves.) 

The  intent  of  the  present  computation  is  to 
show  that  all  three  of  the  isotope  ratios  D/H, 
“C/“C,  and  ”'0/’'0  rnay  well  be  thermometers, 
and  that  their  temperature  coefficients  may  have 
measurable  magnitudes  in  bio-organic  plant  resi¬ 
dues,  and  so  may  give  information  about  climate 
variation  in  the  past. 

For  example,  we  can  measure  these  ratios  in  a 
chronological  sequence  of  tree  lings.  The  results 
of  Parker  [1964]  are  that  (1)  different  plant 
species  growing  side  by  side  in  a  marine  estuary 


TABLU  7.  Temperature  Coefficients  Computed 
for  Isotopes  of  C,  II,  and  O  in  Cellulose 


Fract ionation 

273"  K 

298*k 

Terap. 
Coeff. , 
°/oo/*C 

a|8(°?) 

1.1024 

1 .  18(i5 

♦  0.96 

a  8'COp) 

1 .1200 

1 . 1464 

♦  1.14 

a*H(II^O) 

1.1725 

1 . 1950 

♦0.92 

a'  > 

0.9503 

0.9994 

«0.J0 

a  (IJC09) 

1  . 5881 

i ,sy83 

♦  0.4 

a  (LCOiq 

I.S095 

1.5606 

*2.0 

The  computed  temperature  coefficients  listed 
above  agree  very  well  with  the  value  0.3S  °/oo/*C 
measured  for  13C/I2C  by  Degens  et  at.  (1968),  and 
with  the  value  3  0  °/oo/*C  measured  for  0/11  by 
Schiegl  (I972J. 
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can  have  dtlTrreiit  percentages  of  UC  in  total 
carbon  content,  and  (2)  individuals  of  the  same 
•species  have  about  t he  same  percentage  of  “C. 
Therefore  in  evaluating  past  climates  from  iso¬ 
tope  thermometers,  we  should  measure  changes 
in  a  single  species  of  tree,  and  we  should  measure 
a  chronological  sequence  consisting  of  several 
individuals,  of  the  same  kind  of  tree,  having 
overlapping  life  spans.  By  measuring  trees  of 
overlapping  life  spans,  and  trees  that  grow  in 
sparse  environments,  such  as  mountain  slopes, 
one  may  avoid  spurious  effects  caused  by  the 
relatively  greater  removal  of  CO,  from  air  near 
the  ground  by  abundant  low-growing  plants. 

The  temperature  coefficient  of  oxygen  in  rain 
water  could  be  as  large  as  —0.3(j%c/°C  [Urey, 
1047].  Also  important,  is  the  fact  that  COa 
mixes  rapidly  through  the  entire  global  atmos¬ 
phere  in  a  time  of  ~5  years,  and  with  the  oceans 
[Bien  mul  Suras,  1067;  Stress,  1970;  Craig, 
1057]  in  about  15  years,  so  that  at  any  time 
the  atmospheric  carbon  isotope  ratios  reflect  the 
temperature-dependent  isotope  separation  in  the 
sea,  which  has  coefficients  of  ~0.2I%<,/°C 
for  oxygen  [Bottiuga  and  Craig,  1969],  and 
— 0.11%,/°C  for  carbon  [Emrich  et  al.,  1970], 
One  may  ask  whether  there  is  isotope  ex¬ 
change  between  the  hydrogens  bound  in  old  cel¬ 
lulose  and  in  new  sap.  As  Table  1  shows,  both 
hydrogens  are  bound  very  tightly,  even  more 
so  than  the  C-C  bonds,  so  that,  knowing  that 
radiocarbon  dates  in  tree  cellulose  mainly  agree 
with  toe  tree  ring  dates,  meaning  that  in  heart 
wood  carbon  does  not  exchange,  we  can  hope  the 
same  to  be  true  of  hydrogen  in  heart  wood  But 
this  point  should  be  te.-ted.  When  cellulose  is 
prepared  in  a  thin  layer  or  finely  ground  and 
repeatedly  wetted  and  dried,  there  is  sonic  ex¬ 
change  f Sepal!  anil  Mason,  1961],  but.  those 
conditions  arc  quite  different  from  those  of  cellu¬ 
lose  in  sound  heart  wood. 

Arknou'ledgments.  This  work  was  supported  bv 
the  Advanced  Research  1 ’rejects  Agency  of  the 
Department  of  Defense  and  monitored  by  the  Air 
Korn  Office  of  Scientific  Research  under  contract 
F41620-72-C-0025 

References 

Baxter,  M.  S.,  and  A  Walton,  Fluctuations  of 
atmospheric  C-14  concentrations  during  the  past 
century,  /‘roc.  Hoy.  Sue.  London,  Ser.  A,  Ml, 
105,  1971 

Bien,  G.,  and  II  E.  Suess,  Transfer  and  exchange 


of  C  14  between  the  atmosphere  and  the  surface 
water  of  the  Pacific  Ocean,  paper  presented  at 
the  Symposium  on  Radioactive  Dating  and 
Methods  of  Low  Level  Counting,  U.N  Doc.  SM 
87/55,  Int.  At  Energy  Ag.,  Vienna,  1907 

Bottinga,  Y.,  Calculation  of  fractionation  factors 
for  carbon  and  oxygen  isotopic  exchange  in  the 
system  calcite-earbon  dioxide-water,  J.  Phys 
(’hern.,  72,  800,  1967. 

Bottinga,  Y.,  and  H.  Craig,  Oxygen  isotope  frac¬ 
tionation  between  CO*  and  water  and  the  iso¬ 
topic  composition  of  marine  atmospheric  COi, 
Earth  Planet.  Sri.  Lett.,  S,  303,  1970 

Craig,  II.,  The  Geochemistry  of  the  stable  carbon 
isotopes,  Geochim.  Cosmochim.  Acta,  3,  53 
1953. 

Craig,  H.,  The  natural  distribution  of  radiocarbon 
and  the  exchange  time  of  carbon  dioxide  be¬ 
tween  atmosphere  and  sea,  Tellus,  9,  1,  1957. 

Degens,  E.  T„  R.  II.  L.  Guillard,  W.  M.  Sackrtt, 
and  J.  A.  Hellebust,  Metabolic  fractionation  of 
carbon  isotopes  in  marine  plankton,  1,  Tem¬ 
perature  and  respiration  experiments,  Deep  Sea 
Pea.,  15,  1,  1968. 

Emiliani,  C.,  Paleotemperature  analysis  of  Carib¬ 
bean  sea  cores,  J.  deal.,  74,  109,  1966 

Emrich,  K.,  D.  II.  Ehhalt,  and  J.  C.  Vogel,  Carbon 
isotope  fractionation  during  the  precipitation  of 
calcium  carbonate,  Earth  Planet.  Sci.  Lett.,  8, 
363  1970. 

Glens  an,  J.  D.,  and  I.  Freedman,  Deuterium: 
Natural  variations  used  as  a  biological  tracer, 
Science,  109,  10S5,  1970. 

Herzberg,  G..  Infrared  and  II  man  Spectra,  Van 
Nostrand-Reinhold,  Princeton,  X  J.,  1945. 

Li,  Y.  II..  and  T.  F.  Tsui,  The  solubility  of  COi 
in  water  and  sea  water,  J.  Geophys.  lies.,  76, 
4203.  1971. 

McCrea,  J.  M.,  On  the  isotopic  chemistry  of 
carbonates  and  a  paleotemperature  scule,  J. 
Cheat.  Pl  y-..  1*.  819.  1950 

Merlivat.  L  .  R.  Botter,  and  G.  Nief.  Fraciionne- 
nient  isotopique  au  cours  de  la  distillation  de 
l'eau,  J.  ('him.  Phys.,  CO.  .56.  1963. 

Mook,  \\  .  G..  and  J.  C.  Vogel.  Isotopic  equilib¬ 
rium  between  shells  and  their  environment. 
Science,  150,  S74.  1968. 

Nvdal.  R.,  Variation  tn  Cu  concent  rut  ion  in  the 
atmospherr  during  the  last  several  years,  Tellus, 
I"',  271.  l'.IGG. 

Parker,  R.  L.,  The  bio-geochemi'try  of  the  stable 
isotopes  of  carbon  in  a  marine  bay,  Geochim. 
Cosmochiin.  Acta,  2S,  1155,  1904. 

Rogers,  M.  S„  and  C.  B.  Koons,  Organic  carbon 
delta  Cn  values  from  quaternary  marine  se- 
q lienees  in  the  Gulf  of  Mexico,  a  reflf  tion  of 
palcotem|M  rature  changes,  Trans.  G.dj  Const 
Ass.  Gcrd.  Sue.,  19,  529.  1969 

Sackett,  W  M  ,  \V .  R.  Eckclnianii,  M  1.  Bender, 
and  A.  \V  11.  Be.  Temperature  dependence  of- 
<  arbon  isotope  composition  in  marine  plankton 
and  sediments,  Science,  14$,  235  '965. 

Schiegl,  \V.  E.,  Deuterium  content  peat  as  u 


G-  8 


Multiple  Thermometry 


4317 


paleo-climate  recorder,  Science,  175,  512-513, 
1972 

Sepall,  0.,  and  S.  G.  Mason,  Tritium  exchange 
between  cellulose  and  water,  Can.  J.  Chem.,  39, 
1934,  1961 

Suess,  H.,  Transfer  of  C  14  and  tritium  from  the 


atmosphere  to  the  ocean,  J.  Geophys.  Res..  75, 
2363,  1970. 

Urey,  H.  C.,  Thermodynamic  properties  of  iso¬ 
topic  substances,  J.  Chem.  Soc.,  562,  1947, 

(Received  March  9,  1972.) 


G-9 


